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ABSTRACT 

We report the results of an X-ray observing campaign on the massive, evolved star 77 Car, 
concentrating on the 2003 X-ray minimum as seen by the XMM-Newton observatory. These are 
the first spatially-resolved X-ray monitoring observations of the stellar X-ray spectrum during 
the minimum. The hard X-ray emission, believed to be associated with the collision of r\ Car's 
wind with the wind from a massive companion star, varied strongly in flux on timescales of days, 
but not significantly on timescales of hours. The lowest X-ray flux in the 2 — 1 keV band seen 
by XMM-Newton was only 0.7% of the maximum seen by RXTE just before the X-ray minimum. 
In the latter half of the minimum, the flux increased by a factor of 5 from the lowest observed 
value, indicating that the X-ray minimum has two states. The slope of the X-ray continuum 
above 5 keV did not vary in any observation, which suggests that the electron temperature of 
the hottest plasma associated with the stellar source did not vary significantly at any phase. 
Through the minimum, the absorption to the stellar source increased by a factor of 5 — 10 to Nn 
~3-4xl0 23 cm' 2 . The thermal Fe XXV emission line showed significant excesses on both the 
red and blue sides of the line outside the minimum and exhibited an extreme red excess during 
the minimum. The Fe fluorescence line at 6.4 keV increased in equivalent width from 100 eV 
outside the minimum to 200 eV during the minimum. The small equivalent widths of the Fe 
fluorescence line suggests small fluorescence yield in the companion's low-density wind. The lack 
of variation in the plasma temperature is consistent with the eclipse of the X-ray plasma during 
the minimum, perhaps by a clumpy wind from the primary star, although the deformation of 
the Fe XXV profile and the relatively weak fluorescence Fe line intensity during the minimum 
may suggest an intrinsic fading of the X-ray emissivity. The drop in the colliding wind X-ray 
emission revealed the presence of an additional X-ray component which exhibited no variation on 
timescales of weeks to years. This new component has relatively cool temperature (kT ~1 keV), 
moderate JVh (~5xl0 22 cm -2 ), large intrinsic luminosity (Lx ~T0 34 ergs s _1 ) and a size <1" 
(2300 AU at 2.3 kpc). This component may be produced by the collision of high speed outflows 
at v ~1000— 2000 km s _1 from 77 Car with ambient gas within a few thousand AU from the star. 

Subject headings: Stars: individual (n Car) — stars: early-type — stars: winds, outflows — binaries: 
general — X-rays: stars 
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1. Introduction 

rj Car is a violently unstable, extremely lumi- 
nous object and a key tracer of evolution of stars in 
the up per portion of the Hertzsprung - Russell dia- 
gram ( Davidson fe Humphreys! 1997 ). The star is 
believed to have ha d an initial mass of >150 Mq 
(Hilli er et alJ l200ll ). It is currently in a short, 



poorly understood evolutionary stage, known as 
the Luminous Blue Variable (LBV) phase, which 
is thought to occur near the onset of pulsational 
instabilities. It provides a convenient laboratory 
to study how extremely luminous, massive stars 
evolve and how they shape their environments 
both geometrically and chemically. 

r\ Car is best known for an extraordinarily 
powerful eruption i n 1843 which sent ~12 M Q 
(ISmith et allboOSbh of its atmosphere into space, 
creating a beautiful bipolar nebulosity called 
the Homunculus around the star, rj Car also 
had a minor eruption in 1890, which produced 
a small bipolar nebula inside the Homunculus 



nebu la (the "Little Homunculus" , llshibashi et al 



loss (10" 4 -10- 3 Mp, yr" 1 Davidson et al. 1995: 


Cox et al. 1994 Ishibashi et al.| 


19991 Hillier et al. 


2001; Pittard & Corcoran 12002 


), preferentially in 


the polar direction ( Smith et al. 


2003a: van Boekel et al. 



2003). The ejecta and the stellar wind are rich in 
helium and nitrogen and depleted in oxygen and 
carbon dVerner et al. l2005t iDavidson et alJ 11984 : 



Hillier et al. 200lh. consist ent with CNO process- 



i ng. Davidson et al. ( 2005 ) and Martin fc Koppelmanl 



( 200-d ) suggested from the change of Ha and H/3 



line profiles and the increase of optical brightness 
that rj Car is changing rapidly at the present time. 

The past decade has witnessed an important 
change in our understanding of the star. This 
began with the recognition that the strength 
of some narrow emission lines, notably He I 
10830A, vary predict ably with a period of 5.52 



years (lDamineli| |1996). Al ong with observations 



in infrared (Whitclock et al. 1994; Feast et al 



200 1|) . mm (ICox et al I Il995h. cm dDuncan et al 

' 1999h 



19951 ) , and optical 
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wavelengths, observations of the X-ray emission 
from the star have played a key role. ROSAT ob- 
servations in 1992 first showed a variation which 
appeared to be correlated with Damineli's emis- 
sion line variations (Corcoran et al. 1995), and 
subsequently the 2 — 10 keV lightcurve of the 
star obtained by RXTE (jlshibashi et all Il999h 
showed in detail the X-ray variation culminating 
in a swift, unstable rise to maximum and steep 
fall to minimum lasting for ^3 months, i.e. the 
same time interval as Damineli's spectroscopic 
minimum. These observations clearly show that 
the star varies in a fundamental way every 5.54 
years, and strongly suggest that rj Car is two stars, 
not one. In this binary model, most of the dra- 
matic, pan-chromatic changes are now believed 
to be produced by the interaction of the UV flux 
and wind from a companion star with the wind 
of 77 Car. Variations in the X-ray region are pro- 
duced by a wind-wind collision (WWC) of the 
wind from rj Car with the wind of the companion 
star. A current guess at the system parameters 
describes a massive hot companion with M c ~30 
Mq plus a brighter primary star with M v >80 
M(D in a highly eccentr ic (e > 9), 5.54 year or bit 
( Corcoran et al. I [2001a!, see also Davidson 1999I) . 

During the 1997/98 X-ray minimum, 77 Car 
was the target of every a vailable X-ray obser- 



vatory, including RX TE (llshiba s hi et al 
Corcoran et all l2001af ). ASCA dCorcoran et al 



19991: 



2000), and BeppoSAX (|Viotti et al.l l2002h . The 



ASCA s atellite detected hard X-ray emission from 
rj Car (jCorcoran et al.l 120001 ) . which was char- 
acterized by slightly smaller Ah than the pre- 
minimum state with reduced plasma emission 
measure (E.M.). The limited spatial resolution 
of ASCA, however, left the possibility that the 
observed emission was due to contamination by 
unresolved nearby sources. The BeppoSAX satel- 
lite observed rj Car just after the recovery of the 
X-ray emission ( Viotti et al. 2002! ). The spectrum 
showed strong absorption of Ah~1.5x10 23 cm~ 2 
unlike the ASCA spectrum during the minimum, 
which had an A H ~3xl0 22 cm" 2 . The overall X- 
ray brightness variations observed with RXTE are 
explained well with the colliding wind mechanism 
(|Pittard et all Il998t IPittard fc Corcoran! |2002[) . 
However, key properties of the X-ray lightcurve 
were still poorly understood: dramatic changes 
in X-ray flux as the emission increases to maxi- 
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Table 1 

Logs of the XMM-Newton and Chandra Observations 



Abbr. 


Seq. ID 




Observation Start 
Date 


0ort.it 


Exposure 
(ksec) 


Observation Mode 


Filter 


Before the minimum 


















XMM 00 0726 a 


112580601 


2000 


Jul. 


26, 06:01 (M2) 





470 


• ■ ■ / ■ ■ ■ .29.6 


PFW/PFW/PSW 


thick 


XMM 00 0728 a 


112580701 


2000 


Jul. 


28, 00:51 (M2) 





471 


■■■ / ■■■ /7.8 


PFW/PFW/PSW 


thick 


XMM 03 oi25 b 


145740101 


2003 


Jan. 


25, 13:03 (PN) 





922 


4.8/ ■■■ / ■■■ 


PSW/PFW/PFW 


thick 


XMM 030 127A b 


145740201 


2003 


Jan. 


27, 01:08 (PN) 





922 


4.8/ ■■■ / •■■ 


PSW/PFW/PFW 


thick 


XMM 03 0127B b 


145740301 


2003 


Jan. 


27, 20:42 (PN) 





923 


4.7/ ■•■ / •■• 


PSW/PFW/PFW 


thick 


XMM 030 i29 b 


145740401 


2003 


Jan. 


29, 01:45 (PN) 





923 


5.8/ ■■■ / ■■■ 


PSW/PFW/PFW 


thick 


XMM 03 oi30 b 


145740501 


2003 


Jan. 


30, 00:00 (PN) 





924 


4.8/ ■•■ / •■• 


PSW/PFW/PFW 


thick 


XMM 3060S 


160160101 


2003 


Jun 


08, 13:31 (Ml) 





988 


22.1/27.6/ •■■ 


PSW/PSW/PFW 


thick 


XMM 0306 13 


160160901 


2003 


Jun 


13, 23:52 (Ml) 





990 


21.9/30.4/ •■■ 


PSW/PSW/PFW 


thick 


During the minimum 


















CXOo30720 C 


200216 


2003 


Jul. 


20, 01:47 


1 


008 


90.3 


ACIS-S w/HETG 




XMM 03 0722 C 


145780101 


2003 


Jul. 


22, 01:51 (M2) 


1 


009 


5.8/8.2/8.4 


PSW/PSW/PFW 


thick 


XMM030802 


160560101 


2003 


Aug 


02, 21:01 (M2) 


1 


015 


12.2/17.1/17.5 


PSW/PSW/PFW 


thick 


XMM030809 


160560201 


2003 


Aug 


09, 01:44 (M2) 


1 


018 


8.7/12.3/12.6 


PSW/PSW/PFW 


thick 


XMM030818 


160560301 


2003 


Aug 


18, 15:23 (M2) 


1 


023 


13.0/18.2/18.6 


PSW/PSW/PFW 


thick 


CXOo30828 


200237 


2003 


Aug 


28, 17:38 


1 


028 


18.8 


ACIS-S 




After the minimum 


















CXO030926 


200217 


2003 


Sep. 


26, 22:47 


1 


042 


70.1 


ACIS-S w/HETG 





Note. — Abbr.: abbreviation adopted for each observation. Seq. ID: sequence identification number of each obser- 
vation. Observatio n Start Date: start time of the detector in the parentheses. 4>orbit — (JD [observation start] — 
2450799.792)/2024 ([C or c or an] ||2 00 5) . Exposure: exposure time excluding the detector deadtime; for XMM-Newton ob- 
servations the exposure times arc given for the PN, MOS1 and MOS2 detectors, respectively; The data sets without 
exposure time were not used due to the severe event pile-up. Observation Mode: PSW: prime small window, PFW: 
prime full window. For the XMM-Newton observations the observations modes are given for the PN, MOS1 and MOS2 
detectors, respectively. Filter: optical blocking filter selected for the XMM-Newton detectors. 

a Data from XMM000726 and XMM000728 arc combined and denoted as XMMooo7- 

b Data from XMM030125 > XMM030127A , XMM030127B , XMM030129 and XMM030130 arc combined and denoted as 
XMM0301. 

c Spectra from CXO030720 and XMM030722 are simultaneously fit and those results are denoted as CXOXMMq307- 
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mum, the variation of the absorbing material in 
front of the X-ray emitting region near the X-ray 
minimum, the nature of the rapid fall from X-ray 
maximum to X-ray minimum, the excess in Nn 
after recovery. 

The minimum which occurred in mid-2003 was 
among the best observed astronomical events of 
all time. A key part of this campaign were the 
detailed X-ray observations obtained with RXTE, 
Chandra, and XMM-Newton. RXTE again pro- 
vided crucial monitoring of the daily changes in 
X-ray flux, while Chandra provided previously- 
unobtainable monitoring of the X-ray emission 
line dynamics at key phases of the binary period. 
XMM-Newton obtained critical measurements of 
the spectrum of the source during the low-flux 
state when observations with RXTE or Chandra 
are difficult. 

This paper describes the overall change of the 
X-ray spectrum of rj Car around t he 2003 mini 



1991/92 Decl Janl Febl Marl Aprl Mayl 
1997/96 Jull 



Junl Jull 



Augl Sepl Octl I 
2003 Janl Febl Marl Aprl Mayl 



1 Decl Janl Febl Marl Aprl 

I I I I J 

Junl Jull Augl Sepl Octl 



mum as measured by XMM-Newton (jJansen et al 



20011 ). supplemented wher e necessary with key ob - 
servations with Chandra ( Weisskopf et al. 2.002). 
We try to understand the X-ray emission and ab- 
sorption mechanism and the cause of the X-ray 
minimum, from a comparison of derived spec- 
tral properties such as kT, Ah and emission line 
strengths in these observations and earlier ROSAT 
and ASCA observations during the 1992 and 
1997/98 minima. 

The paper is comprised of the following sec- 
tions. Section 2 describes the observations and 
method of the data reduction, and includes a de- 
scription of problems of unresolved emission com- 
ponents around the central source in the extracted 
XMM-Newton source events. Section 3 describes 
imaging and timing analyses of the XMM-Newton 
data. Section 4 describes the analysis of the X-ray 
spectra during the 2003 minimum. Section 5 com- 
pares the X-ray spectra during the 2003 minimum 
to the previous minima in 1998 and 1992. We dis- 
cuss our results in section 6, and in section 7 we 
summarize our conclusions. 

2. Observations &; Data Reduction 

In 2003, XMM-Newton observed T) Car a total of 
11 times, 5 times in January 2003 (about 6 months 
before the onset of the 2003 minimum), twice in 
June 2003 (weeks before the onset of the 2003 min- 




Orbital Phase 



Fig. 1. — RXTE light curves (solid lines, black: 
2003, red: 1997/98) around the X-ray mini- 
mum with the timings of XMM-Newton, Chandra, 
ASCA, BeppoSAX, and ROSAT pointed observa- 
tions {black 2003, red: 1997/98, green: 1992). 
Multiple ROSAT observations between 1.00 to 
1.01 were summed together for ROSATg 2 o6i2- 



imum, near the X-ray maximum) and four times 
during the 2003 minimum in July and August 
2003 prior to the recovery, which started near 2003 
September 3. Figure Q] displa ys these obser ving 
times on the RXTE lightcurve dCorcoranll2005h for 
context. We also re-analyzed two XMM-Newton 
observations from 2000, whose results were p revi- 
ously published by iLeutenegger et al. ( 2003 ). We 
thus used a total of 13 XMM-Newton observations 
in our analysis, as listed in Table [TJ Throughout 
this paper, individual XMM-Newton observations 
are designated XMM, subscripted with the year, 
month and day of the observation. To improve 
statistics, we combined data from two observations 
in 2000 and combined data from five observations 
in 2003 January because these observations were 
close together in time, and did not show any sig- 
nificant variabilty among them. 

XMM-Newton is composed of three nested 



Wolter I-type X-ray telescopes (jAschenbach et al 



20001 ) with the European Photon Imaging Cam- 
era (EPIC) CCD detectors (PN, MOS1 and 
MOS2) in their fo cal planes (|Striider et all 12001 



Turner et al. I l200ll) . In all the observations, n Car 
was on-axis but the observations were obtained 
with different satellite roll angles. The EPIC 
instrumental modes are listed in Table [T] The 
XMM-Newton spectra suffered significant photon 
pileup in the data taken in prime full window 
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(PFW) mode outside the 2003 minimum when the 
star was bright in the 2 — 10 keV band. We did not 
use these data sets in our analysis. XMM-Newton 
is also equipped with the Reflection Gra ting Spec- 
trometer fRGS; Iden Herder et aill200l[ ). but the 
RGS has limited sensitivity above ~1 keV, where 
emission from the WWC region dominates. We 
therefore did not use the RGS data. 

The analysis of the XMM-Newton data was per- 
formed with version 5.4.1 of the SAS_j] software 
package and version 5.2 of the HEASofil analy- 
sis package. The Observation Data Files (ODF) 
data were processed using the SAS scripts "em- 
chain" and "epchain." We removed events close to 
hot pixels or outside the field of view, and selected 
events with pattern ^4 for EPIC PN spectral anal- 
ysis and pattern ^12 for timing analysis of the 
EPIC PN data and spectral and timing analysis of 
the EPIC MOS data. Most of these observations 
luckily avoided high background periods and the 
instrumental background is negligible in any ob- 
servations within the 0.3—10 keV band. We did 
not need to reject any high background periods 
because the strong emission from r\ Car greatly 
surpassed the observed background levels. 

Due to the limited spatial resolution of XMM- 
Newton, the extracted n Car data inevitably in- 
clude emission from a number of sources (see Fig- 
ure [2]) . The extracted spectrum is usually domi- 
nated by hard X-rays from the wind- wind collision 
in r\ Ca r, which varies from observation to obser - 



vation (jCorcoran et all l200d : IViotti et all [20021. 



The Homunculus Nebula emits weak diffuse re- 
flected X-rays from the central source, which can 
be spatially resolved only with Chandra an d only 
during the minimum (jCorcoran et alj 120041 ). The 
outer ejecta extending about 1' beyond the Ho- 
munculus Nebula emits diffuse non-variable X- 
rays belo w ~1 keV with a ring- or shell-like mor- 
phology (jSeward et allEoOl. IWeis et al.ll2004l ). To 
correct the XMM-Newton spectra for emission 
from the outer ejecta and the Homunculus Nebula, 
we used three Chandra data sets from the 2003 X- 
ray observing campaign, two observations during 
the 2003 minimum and an observation just after 
the 2003 minimum. These Chandra observations 
spatially resolve the outer ejecta emission from the 



emission of the stellar source and so allow us to de- 
termine the spectrum of the outer emission, which 
we can then use to help model the XMM-Newton 
spectra. 

The log of Chandra observations and the tim- 
ing of the observations are also given in Table [1] 
and Figure [T] Chandra observations are desig- 
nated CXO, subscripted with the year, month 
and day of the observation, similar to our des- 
ignation of the XMM-Newton observations. The 
Chandra observations were obtained using the Ad- 
vanced CCD Imaging Spectrometer detector us- 
ing the Spectrometer array (ACIS-S) either with 
or without the high energy transmission grating 
(HETG). For the grating data, Oth-order photon 
events were used for CXO030720 (during the X- 
ray minimum, when the source was too faint at 
1st order), and lst-order photon events were used 
for CXO030926 after the recovery since the central 
source had severe pile-up in the Oth-order data. 
The observation CXO030828 j taken during the min- 
imum with the ACIS-S imaging array with no grat- 
ing, also suffered mild pile- up (~ 4%). The anal- 
ysis of the Chandra data was performed with the 
CIAC| software package, version 2.3 (CALDB ver- 
sion 2.22), version 3.0 (CALDB version 2.23) and 
version 3.0.2 (CALDB version 2.26) for sequences 
200216, 200237, and 200217 respectively. We also 
used version 5.2 of the HEAsoft package, and fol- 
lowed the recommendations of the analysis science 
thread^- A full analysis of the Chandra spectra 
will be published separately (Corcoran et al. 2007, 
in preparation). 

3. X-ray Images and Time Variability 

All XMM-Newton images similarly showed 
emission from the central source and the outer 
ejecta (Figure [2]). Because these XMM-Newton 
images cannot completely separate emission from 
the outer ejecta and central sourc e nor resolve the 



faint X-ray Homunculus Nebula (jCorcoran et al 
2004f ), we used an ellipse of 50" x 37."5 as the 



source region (called "XMM SRC" in Figure ^ , 
which includes all these emission components but 
excludes nearby X-ray luminous stars (such as 
HDE 303308) more than ~ 1' distant. For some 
small window mode MOS observations, the south- 



^-http: / /xmm. vilspa.esa.es/external/xmm_sw_cal/sas_frame.shtml , 
2 http: / /heasarc. gsfc.nasa.gov/docs/software/lheasoft/ 



http:/ /cxc. harvard.edu/ciao/ 
4 http : //cxc . harvard . edu/ciao/threads 



5 



Fig. 2. — True color images of rj Car. left: Chandra HETG Oth-order data from CXO030720 (red: 0.2—1 keV, 
green: 1—3 keV, blue: 3—10 keV); middle: XMM-Newton image accumulated from all the available MOS 
data taken in prime full window mode (red: 0.2—1 keV, green: 1—3 keV, blue: 3—9 keV); right: ASCA GIS 
data from ASCA971224 (red: 0.8-1 keV, green: 1-3 keV, blue: 3-10 keV). The color in the ASCA image 
has been chosen to highlight the diffuse emission. Solid lines show source (SRC) or background (BCD) 
regions (Designations: OE: Outer Ejecta, HN: Homunculus Nebula, CS: Central Source). The Chandra and 
ASCA images were smoothed with a Gaussian with a = 1.5 pixel, while the XMM-Newton image was not 
smoothed. 



east tip of the ellipse fell outside the small win- 
dow. We determined background using regions 
from apparently source-free regions on the same 
CCD chip. 

Table O shows detected photon count rates in 
the soft (0.3-1 keV), medium (1-4 keV), and hard 
(4-10 keV) bands in each detector for each obser- 
vation. The soft band mainly includes emission 
from the outer ejecta, while the medium band in- 
cludes emission from an additional source which 
we call the central constant emission component 
(CCE, see ^4. 1 .2[> . and the medium band out- 
side the minimum and the hard band includes the 
WWC emission (see M4.2|) . In the soft band, the 
MOS2 count rates are 10% lower in XMM 0007 than 
during the 2003 minimum because of the use of 
a smaller MOS2 source region for the XMM0007 
data. The soft-band PN count rates are slightly 
higher outside minimum because the spectral re- 
sponse of the PN allows photons from the bright 
hard source to contaminate the soft band. For 
example, we estimate the soft band contamina- 
tion due to the hard photons in XMM0007 to be 
~0.08 counts s -1 , which makes the corrected PN 
count rates in the soft band 1.06 counts s -1 in this 
observation. This means that the emission from 
the outer ejecta was constant to <2% to better 



than 90% confidence over a half-year interval. On 
the other hand, both the medium and hard band 
count rates decreased dramatically to ~1— 5% of 
the maximum brightness seen by XMM-Newton. 

We constructed light curves with 500 s time 
bins, combining all the available MOS and PN 
data. In Figure 02 we show three sample light 
curves, representing variations before the min- 
imum (XMM030613) and during the minimum 
(XMM030722, XMMososis)- In general, these light 
curves showed no significant variation in most 
energy bands at >90% confidence (Table [2j Fig- 
ure [3]). The exceptions are the medium band light 
curves before the 2003 minimum. The medium 
band light curve from the XMM0306O8 observation 
showed an "excess" of —0.2 counts s -1 (—4 %) 
near the middle of the observation, while the 
medium band light curve from the XMM030613 
observation showed a slight linear increase by 
~ 0.54tg;^ counts s -1 day" 1 (1-cr) or -10% 
day -1 . The linear model is still not accepted 
above the 90% confidence level (x 2 /d.o.f. = 
1.59), due to small fluctuation on timescales of 
<500 s. We note that the hard band light curve in 
XMM030613 is also better reproduced by a linearly- 
increasing model with a slope of — 0.78^ ' 26 
counts s -1 day -1 , or —7.5% day -1 , with x 2 /d.o.f. 
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Fig. 3.— Background subtracted EPIC PN + MOS light curves from XMM 03 o6i3 (top), XMM030722 (bottom 
left), and XMM0308I8 (bottom right). Top and bottom graphs in each panel show the medium (1-4 keV) and 
hard (4-10 keV) bands, respectively. A solid line in each panel shows the best-fit constant model. 
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Table 2 
Time Variability 



Observation Net Count Rate Constant Model Fit 





Soft [0.3-1 keV] 


Med [1-4 keV] 


Hard [4- 10 keV] 




x 2 


/d.o.f. 






PN 


MOS1 


MOS2 


PN 


MOS1 


MOS2 


PN 


MOS1 


MOS2 


Soft 


Med 


Hard 


d.o.f 






(counts s 






counts s 






^counts s 












XMMooot 






0.255" 






1.25" 






0.91" 


0.97 


1.12 


0.82 


76 


XMM0301 


1.15 






7.72 






7.18 






0.98 


1.14 


1.05 


75 


XMM0306O8 


1.10 


0.247 a 




4.58 


1.70" 




6.77 


1.89" 




0.99 


1.45 


0.95 


55 


XMM 03061 3 


1.09 


0.242" 




4.99 


1.84" 




10.49 


2.80" 




1.02 


1.84 


1.29 


62 


XMM030722 


1.03 


0.223" 


0.285 


0.456 


0.170" 


0.188 


0.183 


0.045" 


0.049 


0.81 


0.85 


1.31 


16 


XMM030802 


1.07 


0.222" 


0.288 


0.464 


0.161" 


0.190 


0.517 


0.137" 


0.124 


0.97 


1.21 


0.75 


34 


XMM030S09 


1.06 


0.226" 


0.280 


0.476 


0.178" 


0.188 


0.554 


0.159" 


0.151 


1.13 


0.79 


1.09 


24 


XMM0308I8 


1.05 


0.221" 


0.288 


0.514 


0.169" 


0.202 


0.648 


0.179" 


0.174 


1.34 


0.96 


0.92 


37 



"Source regions arc slightly smaller for these Prime Small Window observations. 



= 0.98. This increase is roughly consistent with 
the RXTElight curve (~15% day" 1 , see Figured]). 



4. X-ray Spectra 

For each XMM-Newton observation, we pro- 
duced EPIC PN and MOS spectra using the same 
source and background regions as in iJ3J Figure 
H presents the EPIC PN spectra for all 2003 ob- 
servations. The X-ray spectra are a combination 
of a variable hard WWC component, along with 
non-variable emission from a number of compo- 
nents. Consistent with Tabled the spectra below 
1 keV are almost constant through the observa- 
tions. Above 1 keV the spectra varied by a factor 
of 2—3 before the 2003 minimum, then decreased 
dramatically by more than an order of magnitude 
after the onset of the 2003 minimum in June, with 
small recovery by a factor of ~3 at E >4 keV dur- 
ing August. Most spectra clearly showed lines of 
hydrogen- and helium- like ions of Mg, Si, S, Ar, 
Ca, Fe, and Ni, and a fluorescence line from cold 
Fe. 

In the following we examine first the non- 
variable emission, then use this emission to fully 
understand changes in the variable component. In 
our analysis, we adop t abundances relativ e to so lar 
abundances given by I Anders fe Grevesse ( 19891 ) 



4.1. Non- Variable Components 

4-. 1.1. Contribution of the Outer Ejecta and Ho- 
munculus Nebula to the X-ray Spectra 

We use Chandra observations to estimate the 
amount of contamination from the outer ejecta 



and the Homunculus Nebula in the XMM-Newton 
spectra. 

We extracted the spectrum of the outer ejecta 
from the Chandra CXO030720 HETG Oth-order 
data, which has the smallest contamination from 
the central hard source. We extracted the source 
spectrum from an 83" x 47" ellipse centered on the 
star ("OE SRC" in Figured]), with its major axis 
parallel to the long axis of the extended outer X- 
ray nebula, but excluded a 20" x 15" ellipse ("HN 
SRC" in Figure^) centered on the central source 
and the Homunculus Nebula. We extracted back- 
ground spectra from a nearby apparently source- 
free region. r\ Car is embedded in soft diffuse emis- 
sion from the Carina nebula, but its contribution 
to the T) Car spectrum is negligible. 

The spectrum from the "OE SRC" region, 
excluding the spectrum from the "HN SRC" 
region, is shown in Figure [5J This spectrum 
shows a strong emission line from nitrogen, which 
re-confirms the N enhancement in the outer 
ejecta found by a number of previou s analy- 
ses, (optical / UV: iDavidson et al.1 1982, X-rays 
Tsuboi et al.lll997l [N/O ~5], iLeutenegger et al 



20031 [N/O >9]1. This spectrum also shows 
a significant excess at 1.3 keV, the energy of 
the helium-like Mg line, which seems to re- 
quire a Mg overabundance of >3 solar at kT 
~0.6 keV (though plasma with that temperature 
should emit a str o nger hydrogen-like Mg line) . 
Leutenegger et al. (I2OO3I) also found a similarly 
strong Mg line in the XMM-Newton grating spec- 
trum, though they conjectured that some of the 
Mg emission might be contamination from the 
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channel energy (keV) 

Fig. 4. — XMM-Newton EPIC PN spectra of r\ Car in 2003. The central source [central constant emission 
(CCE) and colliding wind emission (CWE)] component, the Homunculus Nebula (HN) and the outer ejecta 
components (OE) are marked. Identified emission lines are marked by lines at their rest energies. Fe I marks 
the location of the Fe absorption edge. The red, blue, and grey lines show Chandra best-fit models of the 
outer ejecta (see Figure [SJ , the Homunculus Nebula, and their sum, folded through the EPIC PN effective 
area. 
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stellar X-ray source. Our results suggest that the 
Mg line does in fact originate in the outer neb- 
ulosity. We modeled the extracted, background- 
corrected source spectrum below 3 keV (where 
the emission from the extended source dominates) 
using a simple absorbe d optically thin, thermal 
plasm a model fWABS, Morrison k. McCammonl 
Il983t MEKAL. Mewe et al.l ll995l). including the 
ACTSABS0 component to compensate for the pro- 
gressive low energy degradation of the quantum 
efficiency of the Chandra ACIS detector. We al- 
lowed the abundances to vary, but constrained 
all abundances to have the same value except 
for nitrogen, which was allowed to vary indepen- 
dently. We added a Gaussian with a fixed line 
centroid at 1.31 keV with no intrinsic broaden- 
ing, to account for the helium-like Mg line emis- 
sion. The model successfully reproduces the spec- 
trum of the outer ejecta, yielding a reduced \ 2 
of 1.09 for 81 dof, with the best-fit parameters 
of kT -0.58 keV, AT H -7.2 xlO 20 cm" 2 , with a 
VMEKAL normalization of 7.2xl0~ 4 , elemental 
abundance of —85 solar for nitrogen, near solar 
abundance for the others, and an excess flux of 
—1.5 xl0~ 5 counts cm~ 2 s _1 for Mg. The plasma 
temperature, — 0.6 keV - mainly constrained 
by continuum above 1 keV - agrees well with 
the highest plasma temperature derived from the 
emissi on lines in the XMM-N ewton grating spec- 
trum (jLeutenegger et al. I l2003h . However, the Ah 
is smaller than interstellar ab sorption to r\ Car 
I>3 x 10 



21 cm 2 see §2.2 of iLeutenegger et al 



2003), the nitrogen abundance is overestimated 
from earlier measurements, and the best-fit model 
does not reproduce the Fe L and Ne lines near 
1 keV that are observed with the high resolu- 
tion grating spectrum. This is perhaps caused by 
assuming a simplistic IT model for a spectrum 
which may have multiple cooler plasma compo- 
nents. Therefore, this model reproduces CCD- 
resolution spectra (>100 eV at 1 keV) at energies 
above —0.4 keV but not higher-resolution grating 
spectra. 

The reflected X-ray emission from the Ho- 
munculus Nebula did not vary significantly be- 
tween CXO nan72n and CXOr) 3n828i an interval of 



40 days ( Corcoran et al.ll2004 ) . Since this emission 
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5 http : / /cxc . harvard . edu/cont- 
soft /software / ACIS ABS . 1 . l.html 



0.5 1 2 5 10 

channel energy (keV) 

Fig. 5. — Chandra spectrum of the outer ejecta 
from CXC-030720 (OE SRC - HN SRC in Fig- 
ure [2|). The solid line shows the best-fit absorbed 
IT model and the dashed line depicts the Gaussian 
component for helium-like Mg. 

represents the accumulated X-rays reflected from 
the entire Homunculus Nebula, any time varia- 
tion of the central source shorter than the average 
light-travel time from the star to the reflecting 
site inside the Homunculus Nebula (—88 days) 
will be smeared out. We assume that the emission 
from the Homunculus Nebula did not vary be- 
tween CXO030720 and CXO030828, which includes 
all XMM-Newton observations during the 2003 
minimum (XMM030722 ~ XMMo308is)- To show 
the relative magnitude of the contamination from 
this reflected component, in Figured] we include a 
model fit of the Chandra spectrum of the reflected 
X-rays from the Homunculus Nebula, which we 
extracted from a 20" x 15" ellipse centered on 
the star (HN SRC in Figure [2]), excluding a 2."5 
radius circle (CS SRC). 

Though the X-ray intensity of the Homuncu- 
lus in observations outside the 2003 minimum is 
less certain, at these times the X-ray emission 
from the central source was so bright that the X- 
ray contamination from the Homunculus is <1% 
even if the reflected emission from the Homuncu- 
lus brightened by a factor of —3. 

4-1-2. Discovery of a "Constant" Component 
Near the Central Source 

Figure H depicts the XMM-Newton EPIC PN 
spectra recorded through 2003, along with the 
best-fit spectral models of the outer ejecta, the 
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Table 3 

Best Fit Parameters of the CCE Component 



Parameters 




IT fit 


kT 


(keV) 


1.05 (0.96-1.13) 


Z(Si) 


(solar) 


0.22 (0.14-0.29) 


Z(S) 


(solar) 


0.47 (0.39-0.55) 


E.M. 


(10 56 cm- 3 ) 


7.0 (6.0-8.2) 


N H 


(10 22 cm^ 2 ) 


4.96 (4.77-5.16) 


ix a 


(10 34 ergs s" 1 ) 


1.3 


X 2 /d.o.f 




1.51 


d.o.f 




1205 



Note. — Elemental abundances of the hard com- 
ponent were found to be non-solar, but do not vary 
with observation. 

a Absorption corrected X-ray luminosity between 
0.5—10 keV, assuming the distance of 2.3 kpc. 



X-ray Homunculus Nebula, and their sum. Ex- 
cept for the excess below 1 keV, which is produced 
by the poor absolute flux calibration of the EPIC 
Prnj, the excess over the summed spectrum repre- 
sents emission in a region (CS SRC in Figure [2]) 
within 2. "5 of the hard X-ray peak. Apparent in 
Figure [4] is an excess over the model in the energy 
range 1 — 3 keV. This emission apparently did not 
vary in time with the harder emission, and had 
the same flux in the Chandra spectra of the central 
region in observations CXO030720: CXO030828, and 
CXO030926 (see also Figure^, so that the emission 
was securely constant for more than two months 
between July 20 - September 26 2003. We call this 
component the Central Constant Emission compo- 
nent, or "CCE". 

The CCE shows lines of helium-like S and Si 
ions but no lines of hydrogen-like S and Si ions. 
This indicates that the CCE comes from relatively 
cool plasma. To estimate the emission parameters 
of the CCE, we subtracted contributions of the 
outer ejecta and Homunculus Nebula from both 
the PN and MOS spectra (see the top panel of 
Figure [6] which shows the PN spectra). We ig- 
nored data bins below 1.4 keV, where the MOS 
and PN spectra have significant discrepancies due 
to problems with the relative flux calibrations be- 
tween the two instruments. To better estimate the 
spectrum at low energies, we also fit the Chan- 
dra spectra of the central region from CXO030720 



and CXO030828 (see the bottom panel of Figure^. 
Because of the finer spatial resolution of Chandra, 
these spectra are mostly free from contamination 
by emission from the outer ejecta. Each XMM- 
Newton and Chandra spectrum was fit by a com- 
bination of a low and high temperature, optically 
thin thermal plasma model with independent ab- 
sorbing columns for each component, including a 
Gaussian line to account for the fluorescence iron 
line at 6.4 keV. The low temperature component 
accounts for the emission from the CCE. 

The parameters of the low temperature compo- 
nent were tied among all the spectra, and we as- 
sumed fixed solar abundances for all elements ex- 
cept for silicon and sulfur, whose abundances were 
allowed to vary simultaneously for all the spectra. 
We allowed the parameters of the hard component 
(which represents the variable emission of the stel- 
lar source) to vary among all the spectra, but we 
kept the abundances at the same (non-solar) val- 
ues in each spectrum. Table |3] shows the best-fit 
parameters for the constant component. The spec- 
trum can be fit with kT -1.1 keV, N H -5.0xl0 22 
cm -2 , with an absorption corrected luminosity of 
log Lx —34 ergs s _1 , and low abundances for Si 
(-0.21 solar) and S (-0.47 solar). 

The model, however, was not accepted at above 
90% confidence for the following reasons: i) The 
Chandra spectra show significant excess below 
—1 keV, perhaps originating in foregr ound emis 



see the XMM-Newton calibration report CAL-TN-0018-2-3 



sion from the outer ejecta "bridge" (jWeis et al 
2004 see also Figure [5]) . The excess can be re- 
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niU-Nswtan Spectr& during the X-ray mn<TniiTn 




1 s 
channel energy (luV) 

Fig. 6.— Top: XMM-Newton EPIC 

PN spectra during the 2003 minimum 
(XMM 03 o722-XMMo308i8) after subtracting 
emission from the outer ejecta and Homunculus 
Nebula. Bottom: Chandra spectra of the central 
region in CXO030720 (black, HETG Oth-order) 
and CXO 03 0828 (red, ACIS-S). The solid lines 
show the best-fit model of a simultaneous fitting 
of XMM-Newton (PN and MOS) and Chandra 
spectra by an individually absorbed 2T model. 
Barred, dot-bar, and dotted lines show the CCE 
component, variable component, and a Gaussian 
for the Fe K fluorescence line, respectively. 



produced by including some flux from the outer 
ejecta using the model derived in §4.1.11 This 
slightly improves the reduced \ 2 horn 1.51 for 
1205 to 1.50 for 1203 degrees of freedom. The 
included flux is ^0.4% of the whole outer ejecta 
emission, ii) The model overestimates the spec- 
tra near 2 and 2.7 keV, where lines of hydrogen- 
like Si and S ions are, and the model underes- 
timates the spectra near 1.8 and 2.5 keV, where 
lines of helium-like Si and S ions are. This suggests 
that the constant component includes plasma with 
temperature cooler than kT ~1 keV, though an 
inconsistency near 2 keV between XMM-Newton 
and Chandra spectra could also be caused by cali- 
bration uncertainties of the Chandra effective area 
around the iridium M-edge. Since cooler plasmas 
have lower emissivities in those lines, our simple 
IT model fitting for the CCE component can un- 
derestimate the elemental Si and S abundances, 
iii) The model underestimates the XMM-Newton 
spectra near 3—4 keV, possibly because the vari- 
able component is a multiple temperature plasma 
(see §4.2.2p . iv) The Fe K line profile is complex 
and cannot be fit by a simple 1-temperature equi- 
librium model (see §4.2.3|) . Nevertheless, the pa- 
rameters given in Table [3] should approximate the 
physical properties of the CCE, except perhaps for 
the elemental abundances, which are more depen- 
dent on details of the model. 

As seen in Figure El the emission between 
1—3 keV in CXO030720 comes predominantly from 
the CCE component. The 1—3 keV image in 
CXOo 3 o720 was apparently point-like, which re- 
stricts the projected plasma size to within ~1" 
of rj Car, the width of the Chandra point-spread 
function. This is equivalent to a projected distance 
<2300 AU at d ^2.3 kpc, suggesting that the hot 
gas which produces the CCE component is inside 
the Little Homunculus Nebula (which has a pro- 



jected extent of ~ ±2", llshibashi et alJl2003R 



4.2. The Variable Component 

In order to determine the spectrum of the vari- 
able stellar component, we then subtracted from 
all the original spectra the recognized non- variable 
emission components - the outer ejecta, the X- 
ray Homunculus Nebula, and the CCE compo- 
nent. The remaining emission represents the vari- 
able component associated with the stellar source, 
and is shown in Figure [7] The top panel shows 
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Fig. 7. — Top: XMM-Newton spectra of the variable emission. Bottom: XMM-Newton spectra renormalized 
to the line intensity of the helium-like Fe line. Best-fit models of the variable emission in XMM0007 and 
CXO030926J convolved with the EPIC PN effective area, are also displayed in grey and orange. Each color 
represents the same observation as in Figure |H Boxes in the bottom panel show the high energy slopes of 
the bremsstrahlung emission at kT —A keV with Nn — 0, 2, and 10 xl0 23 cm~ 2 (top), the bremsstrahlung 
emission at kT =1, 2, 4, and 8 keV without absorption (middle), and Fe line profile of the APEC model at 
kT =2, 4, and 8 keV with abundance of 0.5 solar (bottom). 
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the variable component "corrected" for the con- 
stant components we identified. For comparison, 
we overlaid on the top panel of Figure[7]the best-fit 
models for the XMM oo7 and CXO030926 datasets. 

4-2.1. Variation of the Overall Spectral Shape 

The bottom panel of Figure [7] shows the "cor- 
rected" variable components normalized at the 
Fe XXV line energy to highlight changes in the 
shape of the spectra. From XMM0007 to XMM0301 
(an interval of 2.5 years) the X-ray flux increased 
by a factor of 2-3 without remarkable changes 
in spectral shape. Four months later (observa- 
tion XMM0306O8), the X-ray emission decreased 
only below 5 keV. Five days later (observation 
XMM030613), when 77 Car was close to t he peak of 
the la st flare before the 2003 minimum ( Corcoran! 
20051 ). the hard X-rays above ~3 keV increased by 
a factor of two. The X-ray emission dropped by 
almost two orders of magnitude by XMM030722 , 
the first observation after the 2003 minimum. In 
XMM030802 , the hard band flux recovered by a 
factor of 3 without any significant change in the 
soft band. The hard band flux stayed the same 
in XMM030802, XMM 030809 , XMM0308I8, and 
CXO030828 (not shown) , while the soft band emis- 
sion slightly inc reased. The RXTE light curves 
( Corcoranll2005l ) also show a clear transition in the 
flux level and hardness ratio between XMM030722 
and XMM030802 ■ This means that the X-ray mini- 
mum has two states. By CXOo30926, only the hard 
band flux recovered to the pre-minimum level. 

The bottom panel in Figure [7| clearly shows 
that, throughout the observations, the hard band 
slope above 7 keV did not change significantly, nor 
did the ratio of hydrogen-like to helium-like Fe ion 
lines. This means that the electron temperature of 
the hottest plasma did not change during the 2003 
minimum. On the other hand, the soft band flux 
relative to the hard band decreased gradually from 
the X-ray maximum through the minimum to the 
recovery. This looks like an increase of the absorp- 
tion column to the X-ray plasma though the sit- 
uation is somewhat more complicated, as the fol- 
lowing sections describe. The 2003 minimum can 
be better described as an apparent decrease of the 
emission measure ( KM.) as suggested by earlier 
ASCA observations ( Corcoran et al.ll2000n . which 
means that either the amount of X-ray emitting 
material has declined, or that the amount that is 



visible to the observer is smaller. 

4-2.2. 1- Temperature Fit of the Entire Spectrum 

To quantitatively describe the spectral shapes 
and compare them with earlier results obtained by 
ASCA and BeppoSAX (see $5]), we fit each spec- 
trum by an absorbed IT thermal APEC0 model 
(which provided a better fit near the Fe K lines 
than the MEKAL model), with a Gaussian at 
6.4 keV to account for the Fe K fluorescence line. 
The results are listed in Table HI 

Most of the fits were not acceptable at the 
90% confidence level. One reason for the poor 
fit by the absorbed IT models is that the strong 
Fe Ka profile is quite complicated, as described 
in £ 14. 2 .31 Another reason, in particular for the 
spectra before the 2003 minimum, is that the 
plasma is really composed of multiple compo- 
nents with diff e rent te mperatures as shown by 



Corcoran et alj (|2001bf ). as their best-fit plasma 
temperatures kT ~4-5 keV do not account for 
the strong lines of helium-like S and Si ions seen 
in the XMM-Newton spectra (which should be 
emitted from plasma with kT ~ 1 keV). The 
plasma temperatures outside the 2003 minimum 
(kT ~ 4-5 keV) are consistent with earlier re- 



sults (iTsuboi et all Il997t I Corcoran et all Eooot 



IViotti et all 120021 iLeutenegger et all l20Q3h . but 



they are probably overestimated by our IT model 
fit (see the next section). The derived absorption 
column, which does not depend strongly on tem- 
perature for kT ~3— 5 keV, gradually increased 
from Nn « 5 X 10 22 cm~ 2 to Nu « 4 x 10 23 cm~ 2 . 

4.2.3. Fits to the E >5 keV Spectrum 

The EPIC spectra marginally resolved emission 
lines below ~5 keV and could not reliably mea- 
sure the neighboring continuum level. Spectral 
modelling of this energy band with unconstrained 
elemental abundances, therefore, has intrinsic un- 
certainty. On the other hand, the spectrum above 
5 keV is relatively simple and provides an unam- 
biguous measure of the hottest plasma from the 
continuum slope and Fe K line intensities. To 
characterize the hottest plasma, we thus fit the 
spectra above 5 keV by a simple absorbed IT 
model. 



7 http: / /cxc. harvard.edu/atomdb/ 
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Table 4 

The best-fit models for the whole spectra 



Observation 


4>orbit 


kT 


Z 


Flux 6 .4 kcV 




Observed Flux 


reduced x 2 (d.o.f.) 






(kcV) 


(solar) 


(10~ 4 counts cm -2 s" 1 ) 


(10 22 cm" 2 ) 


(10 _11 ergs cm" 2 s" 1 ) 


XMM0007 


0.470-0.471 


4.6 (4.5-4.8) 


0.67 (0.63-0.70) 


0.91 (0.75-1.04) 


5.1 (5.0-5.2) 


6.5 


2.85 (315) 


XMM0301 


0.922-0.924 


4.3 (4.2-4.4) 


0.64 (0.62-0.65) 


2.4 (2.3-2.6) 


4.8 (4.7-4.9) 


15.0 


10.29 (249) 


XMM0306O8 


0.988 


4.5 (4.4-4.6) 


0.83 (0.81-0.85) 


4.1 (3.9-4.2) 


7.4 (7.3-7.5) 


12.8 


12.77 (417) 


XMM030613 


0.990 


5.4 (5.4-5.5) 


0.95 (0.94-0.97) 


6.9 (6.7-7.1) 


9.1 (9.0-9.2) 


19.3 


17.49 (548) 


CXOXMM0307 


1.008-1.009 


3.1 (2.8-3.6) 


0.44 (0.37-0.53) 


0.16 (0.13-0.20) 


28 (25-32) 


0.19 


1.24 (193) 


XMM030802 


1.015 


3.6 (3.3-3.9) 


0.64 (0.56-0.71) 


0.65 (0.58-0.78) b 


42 (39-45) 


0.75 


2.44 (119) 


XMM 30809 


1.018 


2.9 (2.8-3.1) 


0.64 (0.58-0.71) 


0.55 (0.45-0.65) b 


36 (34-38) 


0.82 


1.19 (94) 


XMM0308I8 


1.023 


2.9 (2.8-3.1) 


0.58 (0.54-0.63) 


0.37 (0.26-0.42) b 


29 (27-30) 


0.96 


1.30 (161) 


CXOo30828 a 


1.028 


4.6 (3.4-5.4) 


0.48 (0.37-0.60) 


0.31 (0.3-0.43) 


18 (17-21) 


0.73 


1.54 (137) 


CXOq30926 


1.042 


2.9 (2.3-3.3) 


0.40 (0.36-0.45) 


3.6 (3.0-4.3) 


36 (33-40) 


6.6 


1.40 (242) 



Note. — These fits refer to the spectrum between 2-10 keV after subtraction of the CCE and the Homunculus emission. Parentheses show the range of values 
for which Ax 2 — 2.7 from the best-fit value. 

a About 4% of the photon events arc piled-up, causing the spectrum to appear slightly harder. 

b These values are upper-limits since the line shapes are not clearly seen in the spectra. 
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Fig. 8.— Spectra above 5 keV: EPIC PN (black), EPIC MOS1 (red) and MOS2 (blue), ACIS (green) 
and HETG (grey: +lst order, brown: —1st order). The solid lines show the best-fit absorbed optically 
thin-thermal plasma model (APEC code) with Gaussian components for Fe K a and 0. 
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Fig. 8. — Continued. 



In each observation, we simultaneously fit all 
the available spectra, using an APEC thermal 
equilibrium model for the thermal emission. We 
simultaneously fit the spectra of two observations, 
CXO030720 and XMM030722 , which were obtained 
during the 2003 minimum within 2 days of each 
other. A preliminary fit to the bright phase spec- 
tra showed that most spectra could be fit with 
kT =3.3 keV, a nickel abundance near 0.8 solar, 
with solar abundances for other elements. There- 
fore, we fixed the nickel abundance at 0.8 solar and 
the temperature at 3.3 keV, and we varied the Fe 
abundance, column density and emission measure. 
We included 2 narrow Gaussian lines to account 
for lines of Fe Ka and K/3 fluorescence, fixing 
their centroids at 6.4 keV and 7.1 keV, with their 
line widths (a) fixed at 0.01 keV, and the K/3 line 
flux fixed at 11.3% of the Ka line flux. We used 
ph otoelectric absorption cross-section s calcu lated 
by iBalucinska-Church fe McCammonl (l992). Ta- 
ble [5] gives the best-fit parameters, and the spec- 
tra, along with their best-fit models, are shown 
in Figure [5J The spectra except for XMM oo7 
and CXOXMM0307 accept the assumed model at 



<90% confidence, mostly because these spectra 
show several spectral features which are not fit by 
the model. 

The most prominent of these features are 
excesses on the blue and red sides of the K- 
shell lines of helium-like iron. Those excesses 
are not caused by poor energy or gain cali- 
bration, since there is good consistency in pro- 
file between the XMM030613 spectrum and the 
near-contemporaneous Chandra grating spectrum. 
The blue excess could be explained if the K 
lines of Fe XXV are Doppler broadened, with 
Av ~4000 km s _1 . This is similar to the de- 
rived wind v elocity of the companion star, v w 
3000 km s" 1 (jPittard fe Corcoran|[2002l) . The red 
excess, especially during the 2003 minimum and 
CXOXMM0307 and XMM 03 o802 when the X-ray 
flux was the weakest, would require a Fe XXV 
Doppler shift of At; ^7000 km s _1 , and such high 
velocities are not expected in r\ Car. However, 
the red excess may be produced by unresolved 
emission lines of Fe in ionization stages below 
Fe XXV in the 6.5-6.6 keV band. Interestingly, the 
ASCA minimum spectrum in ASCA971224 showed 
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Table 5 

The best-fit models for the spectra above 5 keV 



Observation 


<t>OTbit 






1 


og E.M. 


EW(Fo) 




reduced x 2 (d.o.f.) 








(solar) 




(cm" 3 


) 


(eV) 


(10 22 cm" 2 ) 




XMM0007 


0.470-0.471 


0.48 


(0.45-0.52) 


58.1 


(58.1- 


-58.2) 


87 (62-99) 


17 (14-20) 


1.43 (86) 


XMM0301 


0.922-0.924 


0.48 


(0.47-0.50) 


58.5 


(58.5- 


-58.5) 


99 (85-104) 


17 (16-19) 


4.07 (82) 


XMM0306O8 


0.988 


0.55 


(0.54-0.57) 


58.5 


(58.5- 


-58.5) 


164 (158-172) 


21 (20-22) 


4.21 (154) 


XMM 030 613 


0.990 


0.44 


(0.43-0.45) 


58.9 


(58.9- 


-58.9) 


119 (114-123) 


36 (35-37) 


3.66 (245) 


CXOXMM0307 


1.008-1.009 


0.38 


(0.30-0.47) 


57.0 


(56.9- 


-57.0) 


-233 (183-325)" 


40 (32-51) 


1.35 (90) 


XMM 30S02 


1.015 


0.50 


(0.46-0.55) 


57.7 


(57.6- 


-57.7) 


<213 a 


53 (48-58) 


1.91 (71) 


XMM030809 


1.018 


0.70 


(0.63-0.77) 


57.5 


(57.4- 


-57.5) 


<257 a 


32 (26-37) 


1.59 (54) 


XMM 03 0818 


1.023 


0.62 


(0.57-0.67) 


57.5 


(57.5- 


-57.5) 


<150 a 


27 (23-30) 


1.65 (84) 


CXOo30828 b 


1.028 


0.33 


(0.24-0.42) 


57.5 


(57.4- 


-57.6) 


~149 (96-214) a 


37 (28-47) 


1.81 (54) 


CXOo30926 


1.042 


0.33 


(0.29-0.37) 


58.7 


(58.6- 


-58.8) 


110 (90-136) 


56 (50-62) 


1.54 (130) 



Note. — These fits refer to the spectrum after subtraction of the CCE and the Homunculus emission. The plasma temperature 
and the Ni abundance were fixed at 3.3 keV and 0.8 solar, respectively. 

"These values are regarded as upper-limits since the spectra did not clearly show the iron fluorescence line peak. 

b About 4% of the photon events pile-up, causing the spectrum slightly hard. 



a similar excess (see N5.2[) . while a Chandra grat- 
ing spectrum near apastron (4> ~0.53) showed a 
similar broad-band excess ( Pittard fc Corcoran 



2002). Presence of lower-ionization iron lines 
would indicate that the hottest plasma is no longer 
in collisional ionization equilibrium at 3.3 keV, 
and that the electron temperature is higher than 
the ion temperature. Non-equilibrium ionization 
(NEI) effects have been claimed in the X-ray spec- 
trum of WR 140, a nother colliding wind system 



(|Pollock et al.ll2005l) 



Similarly, a hump-like spectral feature is per- 
haps present between 5—6.4 keV and becomes 
more noticeable during the early phase of the min- 
ima. There are no specific emission lines in this 
energy band, except a weak Ca XIX line at 4.56 
keV. A similar feature is som etimes seen in the X - 
ray spectra of some AGNs (jTanaka et al.lll995l ). 
produced by gravitationally red-shifted material 
very near the central blackhole. However, there 
is little evidence that rj Car houses a black hole 
because of the lack of short term X-ray variability 
and relatively low X-ray luminosity. 

A second feature is a hard "tail" at >9 keV, 
best seen in the residuals in XMM oo7, XMM0301, 
XMM 03 o608, and XMM 030 6i3 (see Figure H). This 
feature is probably not an instrumental or back- 
ground artifact: none of the spectra (except for 
CXO030828) suffer photon pile-up, and source 
count rates are much higher than the background. 



During the 2003 minimum the tail seems weaker 
or non-existent, though this may simply be an 
artifact since the source is weak at these times. 
The slope above 9 keV in XMM030613 can be fit 
by a bremsstrahlung model w ith kT ^10 keV or 
power-law model with T ~ 2.3. IViotti et al.l ((2004) 



measured a similar photon index for a hard tail ex- 
tending up to 150 keV, seen in a BeppoSAX PDS 
spectrum in June 2000 (<j) ~0.46), but the flux 
between 13—20 keV, 1.4xl0~ n ergs cm~ 2 s _1 , is 
about three times larger than the extrapolation 
of the hard tail we measured in XMMo306i3- The 
hard excess we see in the XMM-Newton data may 
be from an extremely hot plasma in the colliding 
wind region, or from non-thermal emission due 
to lst-order Fermi acceleration at the wind con- 
tact surface, which would produce a population of 
relativistic electrons that can upscatter UV pho- 
tons from the stellar photosphere and/or X-ray 
photons from the wind-wind shock. 

We also clearly see an excess at ~8.5 keV in 
the PN spectra in observations XMM 03 oi and 
XMM 306i3- This feature might be due to emis- 
sion from K shell lines of heavy metals like Cu and 
Zn. Though the EPIC instrument al background 



show s emission from those elements (jStruder et al 
120011 ). inspection of the EPIC background data 
provided by the XMM project showed that the 
instrumental background at this energy was neg- 
ligible in spectra obtained outside the minimum. 
This feature could indicate an overabundance of 
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Cu or Zn, or perhaps be a part of the Cu edge at 
—9 keV (which would also require a large Cu over- 
abundance as well). On the other hand, the fea- 
ture might be produced by absorption edges from 
Fe XXV (- 8.8 keV) and Fe XXVI at -9.3 keV 
(|Lotdll967h . 

Finally, the fourth feature is a marginal dip 
around 7 keV, which is seen in all the PN spec- 
tra, but which is not so prominent in the MOS 
spectra. This dip could be caused by an underesti- 
mate of the line intensity of hydrogen-like Fe ions 
at 6.9 keV (perhaps due to an underestimate of 
the maximum plasma temperature) or an under- 
estimate of the strength of the Fe edge structure at 
7.1 keV, perhaps indicating larger Fe absorption. 
This dip is especially noticeable in the observa- 
tion CXO030926: an d will be discussed more fully 
in the analysis of our Chandra HETGS observa- 
tions (Corcoran et al. in preparation). 

4-2.4- Adding a Soft Component to Reproduce 
the Entire Spectrum 

Extrapolation of the model spectrum above 
5 keV from the XMM 000 7, XMM 030 i, XMM 030 608, 
XMM 03 o6i3, CXC-030828 and CXO030926 data to 
lower energies underestimates the observed emis- 
sion at E <5 keV, indicating the presence of addi- 
tional cooler emitting material. We thus modeled 
the emission below 5 keV by fixing the best-fit 
model spectrum derived in ^4.2.31 and adding an- 
other IT component with independent absorption. 
The results are given in Table [6] and Figure [9] 

These fits were not accepted at >90% confi- 
dence except the fit to the CXO030926 spectra. 
One reason for the poor quality of the fits is that 
the fixed hard component, which reproduces emis- 
sion above 5 keV, somewhat overestimates the 
emission near 5 keV, which might indicate we 
slightly underestimated the absorption to the hard 
component. Another reason is that the best- fit 
models do not reproduce the strengths or loca- 
tions of many emission lines, in particular S XIV, 
which is generally emitted by plasma at kT > 1 
keV. This means that the spectra may require ad- 
ditional components with kT between 1—3 keV. 
The -/Vh values we derive from this 2 component 
modeling are very close to the ones from the IT 
fitting in gTl 



5. Comparison with the Previous Obser- 
vations Near X-ray Minima 

The weekly to daily monitoring of 77 Car with 
RXTE (ICorcoranll2005h showed significant cycle- 
to-cycle variability in flux and X-ray hardness. 
Since the RXTE data are contaminated by in- 
strumental and cosmic background, we compared 
our XMM-Newton and Chandra results with ear- 
lier spatially-resolv ed observati ons from A SCA 
dTanaka et al.lll994h. Be ppoSAX (|Scarsilll993[ ) and 
ROSAT (lTrumperill984h . 

Between phase —0.10 < (f> or Ut ^ 0.10, there 
were 4 ASCA observations and 1 BeppoSAX obser- 
vation near the 1998 minimum, and 1 ROSAT ob- 
servation during the 1992 minimum. These obser- 
vations are summarized in Table [7] Because of the 
extended point responses (FWHM —a few arcmin- 
utes) of ASCA and BeppoSAX compared to XMM- 
Newton, source spectra extracted from these ob- 
servatio ns used relatively large (— 3') extraction 



regions ( Corcoran et al. 2000t Viotti et al. 2002). 



This extraction area is about a factor of —20 larger 
than the area used in the XMM-Newton analy- 
sis (XMM SRC), and about a factor of -5000 
larger than the area used in the Chandra analy- 
sis (CS SRC). Because the neighborhood of r\ Car 
is crowded with X-ray sources and diffuse emis- 
sion, the ASCA and BeppoSAX spectra will suffer 
greater contamination from neighboring sources. 

5.1. Comparison of the 1997/98 and 2003 
Spectra Outside the Minimum 

Outside the minimum, contamination of the 
stellar source spectra in the ASCA and BeppoSAX 
observations from the surrounding sources, the 
CCE component, and the reflection from the Ho- 
munculus Nebula is <10% above —2.5 keV (see 



IViotti et al. ( 2002 ) tried two models to fit their 
BeppoSAX spectra: an absorbed single temper- 
ature bremsstrahlung model + Gaussian for the 
spectrum between 3—10 keV, and an absorbed 
two temperature model to fit the entire spectrum 
between —0.1 — 10 keV. Both gave consistent re- 
sults with each other for the central hard variable 
component. On the other hand, in an analysis of 
ASCA spe ctra obtained ju s t prio r to the 1997/98 



minimum, 



Corcoran et al. ( 2000l ) fixed the emis- 



sion from the outer ejecta using the model derived 
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Fig. 9. — Fits of the spectra below 5 keV. 
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Table 6 

Additional soft component for the spectrum below 5 keV 



Observation 


<t>orbit 




kT 






Z 




log E.M. 




N H 




Reduced x 2 (d.o.f.) 








(keV) 




(solar) 




(cm" 3 


) 


(10 


22 cm" 2 ) 




XMM 000 7 


0.470-0.471 


1.1 


(1.1- 


1.2) 


0.25 


(0.20- 


0.31) 


58.3 (58.2- 


58.4) 


5.4 


(5.2- 


-5.5) 


2.04 (225) 


XMM0301 


0.922-0.924 


1.1 


(1.1- 


1.1) 


0.18 


(0.16- 


0.19) 


58.8 (58.8- 


58.8) 


5.4 


(5.3- 


-5.4) 


8.85 (164) 


XMM 030 608 


0.988 


1.1 


(1.1- 


1.1) 


0.24 


(0.22- 


0.26) 


58.6 (58.6- 


58.7) 


6.9 


(6.8- 


-7.0) 


5.41 (260) 


XMM030613 


0.990 


1.6 


(1.5- 


1.6) 


0.25 


(0.22- 


0.28) 


58.6 (58.5- 


58.6) 


7.7 


(7.5- 


-7.8) 


12.02 (300) 


CXOo30S2S 


1.028 


0.5 


(0.4- 


0.9) 




<0.58 




59.9 (58.1- 


60.1) 


29 


(17- 


39) 


1.32 (80) 


CXO030926 


1.042 


0.9 


(0.6- 


1.8) 


0.28 


(0.05- 


0.75) 


59.2 (58.0- 


60.5) 


23 


(13- 


33) 


0.52 (109) 



Note. — These fits refer to the spectrum after subtraction of the CCE and the Homunculus emission. 
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Orbital Phase 



Orbital Phase 



Fig. 10. — Variation of the temperature and absorption column (TVh [whole]) vs. orbital phase. The filled 
marks are from the 2003 cycle (circle: XMM-Newton, triangle: Chandra, square: simultaneous fit between 
XMM-Newton and Chandra). The open marks are from the 1997/98 cycle (circle: ASCA, triangle: Bep- 
poSAX). The ASCA spectrum during the 2003 minimum at phase 1.004 is consistent with the JVh and kT 
values from CXOXMM0307, as shown. 



Table 7 

Observations around the Minimum From Previous Cycles 



Abbreviation 


Satellite 


Date 


PorMt 


Exposure 




kT 


Rcferecncc 












(kscc) 


(10 22 cm" 2 ) 


(keV) 




ROSAT920612 


ROSAT 


1992 Jun. 


12, 22:33 


-0.998 


23.6 








ASCA970703 


ASCA 


1997 Jul 


3, 09:58 


-0.082 


12.6 


4.16±0.23 


3.79±0.31 


1 


ASCA970719 


ASCA 


1997 Jul. 


19, 15:11 


-0.074 


12.8 


4.09±0.30 


4.12±0.39 


1 


ASCA971224 


ASCA 


1997 Dec. 


24, 09:29 


0.004 


58.8 






1 


SAX9 8 0318 


BeppoSAX 


1998 Mar. 


18, 05:01 


0.045 


39. a 


15.4±0.4 


4.35±0.15 


2 


ASCA9807I6 


ASCA 


1998 Jul. 


16, 04:40 


0.104 


13.4 


5.00±0.52 


4.89±0.88 


1 



NOTE. — JVh and kT refer to the h ottest compone n ts in 2-t emperature f i ts of t he entire spectra; ASCA results refer to 
analysis of the GIS2 spectra. Ref. HCorcoran et al.l H2000D , 2 IViotti et al.l (1200211 

a Exposure time for the MECS23 detector. 



in ICorcoran et al.l ( 1998 ). Though the coolest 
temperature they used to fit the ASCA spectra 
(kT ~0.3 keV) is a bit lower than our adopted 
temperature (kT ~ 0.6 keV) and the resulting dif- 
ference in the outer ejecta flux is about a factor of 
three at ~2 keV, this only produces a ~1% uncer- 
tainty in the flux of the variable emission outside 
the 2003 minimum. The results of their BeppoSAX 
and ASCA analyses, therefore, are comparable to 
the result we derived from our IT fits to the entire 
XMM-Newton and Chandra spectra in §4.2.21 

The plasma temperatures in the 1997/98 ob- 
servations were always around ^4—5 keV, which 
is similar to our results. The derived values of Nn 



varied significantly with phase (Figure [TU)l . reach- 
ing a maximum during the X-ray minimum, and 
declining thereafter. The variation in 2003 was 
similar to the Nr variation in 1997/98, except 
for the interval just at the end of the minimum. 
Though BeppoSAX and Chandra observed rj Car 
at a similar orbital phase (0sAX98O3i8 ~0.045, 
0CXOO3O926 ~ 1-042), the TVh measured with Chan- 
dra in 200 3 is a factor of two larger than the value 
derived bv lViotti et~all (|2002l) from the BeppoSAX 
observations in 1998. This is consistent with the 
behavior of the hardness ratio after the recovery 
as measured by RXTE, which was hig her in the 
2003 recovery than in the 1998 recovery (jCorcoran 
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- 1 2 5 10 - 0.5 1 2 

channel energy (keV) channel energy (keV) 

Fig. 11.— Left. ASCA GIS2+3 spectrum on 1997 Dec. 24. Right: ROSATPSPC spectrum on 1992 Jun 12. 
In each plot the spectrum of r\ Car during the X-ray minimum in 2003 is shown (red: the outer ejecta, blue: 
X-ray Homunculus Nebula, green: CCE component, purple: variable emission in CXOXMM0307; orange: 
surrounding sources within 3'. The solid black line shows the sum of these emission components. 



2005) 



5.2. 



Comparison of the 1997/98 and 2003 
Spectra During the Minima 



In order to compare this dataset with the X- 
ray observations during the 2003 minimum, we re- 
analyzed the ASCA971224 CIS data. We screened 
the "revision2" data taken from the HEASARC 
archival using gisclean and filtered it with the 
standard criteria, which excludes data affected 
by the South Atlantic Anomaly, Earth occulta- 
tion and high background in regions with low ge- 
omagnetic rigidities. For both GIS2 and GIS3 
detectors, we used a 3' radius circle centered on 
r] Car as the source region "3' SRC" , and used as 
background emission from a 3' radius circle des- 
ignated "3' BGD" centered to the north-east of 
rj Car (Figure[2]) where the soft diffuse X-ray emis- 
sion from the Carina N ebula is apparently weak 
( Hamaguchi et all 20061 ). We used the standard 
GIS response files, version 4.0 (g[23]v4_0.rmf) and 
generated ancillary response functions with as- 
caarfver. 3.10. We then merged both GIS2 and 
GIS3 spectra together to improve signal to noise. 
The left panel of Figure [Tl] shows the GIS spec- 
trum from ASCAg7i224- 

Because the stellar source was faint during the 
minimum and the ASCA point spread function 



3 http: / /heasarc. gsfc.nasa.gov/W3Browse/ 



is large, contamination by nearby X-ray sources 
is significant for this observation. We therefore 
tried to estimate the contamination from emis- 
sion from sources within 3' of r] Car "3' SRC" us- 
ing XMM030722 (when the central point source of 
77 Car was the weakest) excluding the source re- 
gion "XMM SRC" used for the analysis of 77 Car 
from the XMM-Newton MOS2 data. We used the 
same background region as we used to estimate the 
ASCA background emission, "3' BGD". About 
10% of the emission from rj Car (which amounts 
to ^20% of the emission from the surrounding 
sources even in XMM030722) contaminates the "3' 
SRC" region. We simulated the contamination 
due to rj Car and subtracted it from the ASCA 
source spectrum. The subtracted spectrum (Fig- 
ure [T2")) does not show any evidence of the ex- 
tremely strong nitrogen line at ^0.5 keV from the 
outer ejecta, suggesting that contamination is well 
removed. 

The resulting spectrum, after subtracting the 
simulated rj Car spectrum, shows weak emis- 
sion lines at ~0.9 and 1 keV from Ne IX and 
Ne X, ~0.8 keV possibly from Fe, and a marginal 
excess around ~3-4 keV from Ca. This spec- 
trum can be fit either by an optically thin, ther- 
mal (MEKAL) model with 2-temperature com- 
ponents or by a cool component plus a power- 
law component for the hard emission at slightly 
below the 90% confidence, assuming common 
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Table 8 

The spectral fit of the surrounding region 





Soft Component 




Hard Component 








Model 


kT log EM. 


kT 


log EM. r 


Z 




reduced x 2 (d.o.f.) 




(keV) (cm -3 ) 


(keV) 


(cm- 3 ) 


(10 -1 solar) 


(10 21 cm" 2 ) 




2T 


0.24 (0.18-0.29) 57.7 (57.2-58.5) 


2.6 (1.9-3.3) 


56.2 (56.2-56.4) 


0.36 (0.19-0. 70) a 


3.6 (2.6-5.5) 


1.28 (92) 


IT + pi 


0.27 (0.23-0.30) 56.5 (56.3-56.8) 




2.5 (2.4-2.7) 


3 (fix) b 


2.6 (2.1-3.3) 


1.26 (93) 



a The abundace parameters of both soft and hard components are tied together. 
b The absorption is insensitive to abundance in the range between and 10 solar. 



Note. — Values in parentheses denote 90% confidence limits. 




0.5 1 2 5 10 

channel energy (keV) 

Fig. 12.— XMM-Newton MOS2 spectrum from 
XMM030722 extracted from a 3' radius circle ex- 
cluding emission from r\ Car. The solid line shows 
the best-fit 2T MEKAL model with common ab- 
sorption, while the broken lines show the individ- 
ual fit components. 



absorption for both components (Table [8]). In 
either model, the average absorption to these 
sources is consistent with the interstellar value 
f~3x 10 21 cm" 2 



see §2.2 of Leutenegger et al 



2003). As the spectrum does not have strong 



emission lines, the elemental abundance of the 2T 
model is small (^0.04 solar), while the abundances 
are not constrained by the IT plus power-law 
model. We thus fixed the abundances at 0.3 so- 
lar, a typical value for stellar X-ray emission ob- 
tained from low resolution CCD s pectra (e.g. OB 



stars , iKitamoto fc Mukai 1996 : Kitamoto et al 
20001: low-m ass MSs. iTagliaferri et al.lll997t low 



mass PMSs. | Yamauchi et al] 1996 : Kamata et al 



19971 : iTsuboi et al.lll998l 1 9 l The observed flux is 

1 (0.5-10 keV). 
includes at least 55 



^3.2x10 12 ergs cm 

The "3' SRC" region 
Chandra X-ray point sources (|Evans et al.ll2003h : 
6 OB stars and 49 unidentified, weak sources 
that might be low mass pre-main-sequence stars. 
The combined absorption corrected flux from 
these source s calcu lated from Table 3 and 4 in 
lEvans et all ([20031 ) is 1.3xl0~ 12 ergs cm" 2 s" 1 
(0.5—2.04 keV). There is also diffuse emission 



9 High resolution grating spectra with Chandra and XMM- 
Newton showed complex line profiles especially from OB 
stars, and therefore the abundance obtained from low res- 
olutio n spectra may not reflect real elementa l abundance, 
(e.g ICassinelli et al.ll200ll ; iKahn et al.ll200lHMiller et al.l 
2002) 



within the source region wh ich is also calcu- 
(|2003l ) to have a flux of 
The total absorp- 
, is 



ergs cm 2 s 



lated from lEvans et al 
~5.7xl0" 12 ergs cm~ 2 s 
tion corrected flux, ~7xl0~ 12 ergs cm" 2 s 
roughly consistent with our result 8.2x 10" 12 
(0.5—2.04 keV, unabsorbed). Moreover, the com- 
bined absorption corrected fluxes of the point 
sources are as large as the flux from the hard 
component (1.7xl0 -12 ergs cm" 2 s" 1 ), and the 
flux from the diffuse emission is as large as the 
soft component, 6.5xl0 -12 ergs cm" 2 s _1 , in the 
0.5 — 2.0 keV band. As seen in Figure [21 many 
point sources are detectable at intermediate en- 
ergies, where the hard component is dominant, 
but they are not clear in the soft band (except 
for HDE 303308). The hard component proba- 
bly represents emission from point sources with 
the soft component dominated by diffuse emission 
and HDE 303308. 

We have compared the ASCA spectrum to a 
model derived from our fits to the Chandra and 
XMM-Newton spectra including emission from 
the outer ejecta, from the Homunculus Neb- 
ula, from the CCE component and the variable 
model which fits the spectrum obtained from the 
CXOXMM0307 observation, along with emission 
from the surrounding sources (diffuse emission and 
point sources) which fall within the 3' extraction 
circle. The ASCA spectrum agrees well with this 
combined model. Assuming that emission from 
the outer ejecta, X-ray Homunculus Nebula and 
the surrounding sources did not vary between 1997 
and 2003 and that spectral shapes of the CCE and 
variable components did not change dramatically 
the flux of both CCE and variable components are 
the same within 50%. 

5.3. Comparison with the 1992 Minimum 

There is one ROSAT PSPC observation of 



rt Car during the minimum in mid-1992. ICorcoran et al 
(|1995l ) used this observation to show that the X- 
ray emission in the hard ROSATb&nd (E >1.6 keV) 
decreased by a factor of 2 at that time. ROSAT 
however had very little effective area above 2 keV 
so that the bulk of the variable emission was not 
observable by ROS AT. 

( 19951) 



Corcoran et al 



extracted a spectrum 
from a 1/85 radius circle excluding a 30" radius 
circle around the nearby bright star HDE 303308. 
To minimize contamination from surrounding 
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sources, we re-extracted the spectrum using an 
elliptical source region of 65" x 47. "5, which 
includes emission from the outer ejecta but 
which excludes HDE 303308 and point and 
diffuse sources around n Car. We took the 
same background region as used in our anal- 
ysis of ASCAg7i224- In analyzing the PSPC 
spectrum, we used the standard response file 
pspcb_gain2_256.rmf from the ROSAT calibra- 
tion databbase, but we generated the ancillary 
response file using the pcarf v. 2.1.3 as appropri- 
ate for this observation. 

The PSPC spectrum is shown in the right panel 
of Figure [TTJ We have overlaid all the spectral 
components of 77 Car on this spectrum as we did 
for the ASCA971224 spectrum, except for the com- 
ponent due to the surrounding sources, whose con- 
tribution would be small in the ROSAT spectrum 
due to the smaller PSPC extraction region. Emis- 
sion between 0.5—1 keV, where the outer ejecta 
component is dominant, is about a factor of 30% 
larger than the model. This difference may be 
due to calibration uncertainties between ROSAT 
and XMM-Newton, or it might be produced by 
another variable source within the ROSAT extrac- 
tion area. Though the XMM-Newton model of the 
CCE component only contributes to the few high- 
est channels of the PSPC detector, this component 
is consistent with the ROSAT spectrum. Thus it 
appears that ROSAT did not detect the stellar 
component at all during the 1992 minimum. 

6. Discussion 

6.1. Emission from the Variable Compo- 
nent 

RXTE clearly showed that the variable emis- 
sion is periodic and leaves little doubt that the 
emission is produced by the WWC in a binary 
system ( Corcoranl 20051 ) . and if so, that the X- 
ray minimum occurs near periastron passage. Our 
results further strengthen this conclusion. The 
temperature of the hottest plasma, ^3 keV, was 
mostly unchanged through the cycle, consistent 
with the collision of stellar winds at te rminal ve- 
locities (e.g., see Ishibashi et al.l fl999h . Further- 



more, the spectral shape varied in a periodic way, 
as well. However, the JVh measured just after the 
end of the X-ray minimum was about a factor of 
two larger during the 2003 cycle than it was during 



the 1998 cycle. 

Most of the light curves derived from the Chan- 
dra and XMM-Newton data showed no short-term 
variations, suggesting that the X-ray plasma is 
produced by a steady source and is probably larger 



than V s ho 



x t P 



where V s hock is the pre-shock 



velocity and t exp is the exposure time of the obser- 
vation. For V s hock ~ 3000 km s" 1 (corresponding 
to shock temperatures of ~ 3 — 4 keV) and expo- 
sure times of ~ 10 ksec (corresponding roughly to 
the XMM-Newton observing time) , the size of the 
emitting region is probably > 0.2 AU, and per- 
haps considerably larger. The RXTE lightcurve 
also shows clear variations in observations sepa- 
rated by as little as 1 day, suggesting that the size 
of the emitting region is 0.2 < r < 2 AU. Although 
the orbital elements are currently uncertain, the 
stellar se paration at periastron may be as little as 



1.5 AU ([Corcoran et alJ l2001aT >. so that the size 



of the emitting region may be comparable to the 
stellar separation at periastron passage. 

Spectra outside the 2003 minimum require at 
least 2 temperatures, kT ^3 keV and 1-1.5 keV, 
consistent wi th the analysis o f a Cha ndra grating 
spectrum by ICorcoran et al. ( 2001b ). The fail- 
ure to fit most of the spectra by 2 temperature 
models in £14.2.41 is perhaps an indication of the 
presence of more than 2 temperatures. This is 
consistent with a hydrodynamic simula tion of the 
WWC plasma (|Pittard fc Stevendll997h . in which 
the plasma temperature is the highest near the 
stagnation point where the winds collide head-on, 
and cooler farther along the bow-shock. The ex- 
cess at the blue side of the Fe XXV line discussed 



in §4.2.31 may represent streams from the stagna- 
tion point flowing outward near the companion's 
wind terminal velocity. 

The RXTE light curves also displayed quasi- 
periodic X-ray "fla res" , which correspond to hard- 
ness r atio maxima ( Corcoran et al.fl997l : ICorcoran 
20051 ). XMM0306O8 occurred near the bottom of 
one of these flares, while XMM030613 occurred near 
a flare peak. A comparison of these spectra in 
Figure [7] shows that the hard band flux went up 
by a factor of two, while the soft band flux be- 
low ~3 keV did not change except that emission 
lines are apparently weaker in XMMo306i3- This 
is consistent with an increase of the hardness ra- 
tio measured with RXTE. This indicates either a 
slight increase in kT and/or E.M. for the hottest 
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plasma, possibly caused by increase of the density 
of the high-temperature gas near the stagnation 
point, as discussed for a similar e vent in the pre- 



vious cycle observed with ASCA (jCorcoran et al 
2000). The emission at E < 3 keV, on the other 



hand, did not vary strongly between those obser- 
vations. One possibility is that hard X-ray emis- 
sion from the highest-temperature plasma at the 
stagnation point of the shock cone varied strongly 
due to density fluctuations of the primary or sec- 
ondary wind, with the soft X-ray emission coming 
from a larger area which averages out such fluc- 
tuations. Another possibility is that the flares are 
related to instabilities in the WWC plasma. 

6.2. The Absorber around the WWC 
Plasma 

Before the 2003 minimum, ATjj measured from 
our fits to the spectra above 5 keV (hereafter 
Ah[>5 keV]) was always about a factor of four 
larger than Nn measured from the spectra be- 
low 5 keV (A H [<5 keV]; Figure [13]). A H [>5 keV] 
mainly represents Fe absorption, while A hJ<5 keV] 



repre sents absorption by lighter elements (jMorrison fe Mc 
Il983l ). Though this discrepancy can be produced 
by an Fe overabundance in the 77 Car ejecta, the 
elemental abundances around the stellar source 
are near solar except for He and N, which are 
overabun dant, and O and C , which are strongly 



that extra absorbing material must be located 
near the hottest plasma. 

After XMM0301 and through the onset of the 
2003 minimum, N-r[>5 keV] and the absorption 
column derived from fitting the spectra over the 
entire 1 — 10 keV band (A^fwhole], see Figure [TO]) 
increased to 3— 5xl0 23 cm' 2 . As shown in Fig- 
ure [10] the observed maximum column density oc- 
curred near the mid-point of the X-ray minimum, 
and the column densities declined thereafter to- 
wards the end of the minimum. During the X-ray 
minimum, Nn[<5 keV] was unmeasurable because 
of the near total suppression of the low energy flux. 
This could mean that the E.M. of the cool compo- 
nent decreased during the minimum, or that the 
absorption to the cool plasma also increased on 
the line of sight. Nevertheless, the A^h[>5 keV] 
confirms that the X-ray minimum is closely asso- 
ciated with an increase in the amount of absorp- 
tion along the line of sight. This suggests that the 
WWC plasma entered into, or was hidden behind, 
the densest part of wind from the primary star. 
The behavior of the absorption towards the hard 
:nt, #h[>5 keV], is in interesting contrast 



depleted (iHillier et all l200ll: IVerner et all [2005; 



Davidson et al.l 1984) . Abundance anomalies of 



He, N, O and C would not change the amount of 
absorption above ~1 keV. This means that the 
soft emission does not suffer as strong absorp- 
tion as the hard emission before the minimum, 
and hence this supports the idea that the vari- 
able emission originates in multiple-temperature 
plasmas with various absorptions. This is con- 
sistent with the idea that the hot emission origi- 
nates from a highly absorbed region near the apex 
of the shock cone, with the lower-energy emis- 
sion arising from a region farther along the shock 
interface which suffers less absorption. Neither 
-/Vh[>5 keV] nor A^h[<5 keV] varied significantly 
between XMM0007 at <j) ~0.47 (close to apastron) 
and XMM0301 at cj> ~0.92, when the stellar separa- 
tion changed by about a factor of 3. The value of 
iVn[<5 keV] before the minimum (~5xl0 22 cm~ 2 ) 
may be absorption beyond the wind of the pri- 
mary. The higher value of JVh[>5 keV] suggests 



to the behavior of Nn [whole]. A^h[>5 keV] had al- 
ready increased to ~4xl0 23 cm -2 by XMM030613 
(i.e. prior to the start of the X-ray minimum) and 
it did not change strongly after that. This could 
be explained if the WWC region was bent, perhaps 
by the Coriolis force, causing the stagnation point 
to enter into the densest part of the primary's wind 
earlier than the downstream material. 

The column density did not decline until some 
time after the 2003 minimum ended. This con- 
firms that the column density variation is not 
completely synchronous with the X-ray minimum. 
Conjunction probably occurred during the interval 
when JVh was near its maximum observed value, 
in the phase interval 0.99 < <fi < 1.05. We note 
that the Nr measured in 2003 (0=1.042) is a fac- 
tor of two l arger than the Nn me asured in 1998 
(0=0.045). iDavidson et~all (|2005l see their Fig- 
ure 2c) also noted that the absorption component 
in the P-r;(A9017) line in 2003.72 (0 w 1.0) was 
stronger than in 1998.21 (<j> « 0.0), which might 
also indicate that the amount of absorbing mate- 
rial in the line of sight near the start of the X-ray 
minimum is increasing with time. However they 
also note that the equivalent widths of the H-a 
emission lines were lower just after the end of the 
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Fig. 13. — Variation of JVh, E.M., intensity and EW of the Fe fluorescence line and RXTE flux in the 
2003 cycle. The XMM-Newton results are shown as circles, Chandra results as triangles, while squares show 
results derived from a simultaneous fit between XMM-Newton and Chandra. The dash line in the top panel 
represents our estimate of the absorption beyond the central source. 
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2003 minimum. This might suggest a distribution 
of wind material which has preferentially grown 
denser in the line of sight than perpendicular to 
the line of sight in the 1998—2003 interval. 

The equivalent width (EW) of the Fe K fluo- 
rescence line varied between 100—200 eV, which 
is far smaller than the EW (~1.5 keV) produced 
by s cattering X-rays from the Homunculus Neb- 



ula ([Corcoran et al.ll2004f ). This suggests that the 



X-ray emission in all the observations comes di- 
rectly from the WWC plasma, even during the 
X-ray minimum, with only a s mall sca t tered com- 
ponent (see Figure E]). From Ihiouel (1985), the 
Fe K fluorescence line EW produced by scattering 
from a spherically-symmetric medium around an 
X-ray source is EW -7.5x7V H [>5 keV]/10 22 eV, 
where we have multiplied Inoue's coefficient by 
0.75, the flux ratio in the 7.1 — 9 keV band be- 
tween a kT ^3.3 keV thermal emission component 
appropriate to the rj Car spectra and the T = 1.1 
power -law emission component assumed by llnoue 
([19851) . In Figure [131 we plot the expected Fe- 
K fluorescence line EWs as open symbols for a 
spherically-symmetric scattering medium. The 
differences between the expected and observed 
equivalent widths seem to be larger when the X- 
ray flux goes up and smaller when it goes down 
(XMM 03 o608, XMM 03 o6i3 and CXO 03 0926), which 
may suggest that the variation of the fluorescence 
line is delayed by light travel time. 

The observed EWs are ^25% smaller than 
the expected EWs even during the XMM0007 
(0=0.470) and XMM 030 i (0=0.923) observations 
when the direct X-ray emission did not vary 
strongly. This suggests that the scattering ma- 
terial i s not spherically symmetric as assumed in 
Inouel (1985), and that there is a deficit of scat- 
tering material off the line of sight. This deficit 
is probably produced by the wind-wind collision, 
since the region on the companion side of the 
shock has a much lower density than on the pri- 
mary star's side. In order to reduce the fluorescent 
Fe emission by 25%, the full opening angle of the 
wind shock cone would need to be ~120°. This 
is consistent with the meas urement of the half- 
openin g angle of 50 — 120° bv lPittard fe Corcoran 
(|2002h . 



6.3. The Nature of the X-ray Minimum 

The 2003 X-ray minimum can be described 
as an apparent decrease in E.M., along with an 
increase in TVfj. However, a column density of 
TVh ^5x10 24 cm -2 is needed to reduce the emis- 
sion at 10 keV by two orders of magnitudes, 
and the observed maximum column density (Na 
~ 5 x 10 23 cm~ 2 ) is much lower than this. The 
ASCA observations du ring the 1997/98 minimum 
( Corcoran et al ] l2000h also appeared consistent 
with an apparent decline in E.M. 

Two mechanisms could produce the observed 
decrease in E.M. One possibility is an "eclipse 
model" , namely that emitting region is almost to- 
tally, but not completely, obscured by an opti- 
cally thick absorber with Nn> 10 25 cm~ 2 , and 
we only see the least absorbed residual emi ssion 
during the min i mum. Ishibashi et al. (1999) and 



Corcoran et al.1 (|2001aD tested the eclipse model 
using the RXTE light curve of rj Car and showed 
that the duration of the X-ray minimum (~ 3 
months) at high energies (E > 5 keV) was difficult 
to explain for a spherically symmetric wind from 
r) Car if the mass loss rate was constant at ~ 10~ 4 
M yr _1 . The eclipse model seems to require that 
the wind is not spherically symmetric, or that the 
absorption to the hard X-ray emitting region was 
enhanced for an extended interval near periastron 
passage. It also suggests that a portion of the hard 
emission is visible during the minimum so that the 
size of the hottest region is greater than the size 
of the occulting region. 

Another possibility is that the emissivity of the 
WWC itself fades during the minimum, i.e. that 
there is a decline in the amount of material hot 
enough to generate ~ 5 keV X-rays. Some WWC 
models predict a change in the X-ray emissiv- 
ity of the shock around periastron. For ex am- 
ple, "radiative braking" (jGavlev et al.l Il997l ) in 
which UV photons from r\ Car decelerate the com- 
panion's wind at phases near periastron could re- 
duce the intrinsic X-ray emissio n from the wind- 
wind collision at E > 5 keV. iDavidsonl (|2002h 
suggested that strong instabilities near periastron 
might cause the shocked gas to radiat ively cool, re - 
ducin g the emissio n at hi gh energies. ISoker ([20051) 
and Akashi et al. ([20061) suggested that Bondi- 
Hoyle accretion of the primary star's wind by the 
companion star near periastron could suppress the 
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Fig. 14. — Cartoon of the r\ Car system. Winds from the primary star and the companion collide together, 
to form the wind-wind collision (WWC) region which has a conical shape. Behind the WWC boundary is a 
low density region filled with the wind of the companion star. The apex where the winds collide head on has 
hotter gas than the downstream material and also accumulates a buildup of cool primary wind material in 
front. We view the X-ray emission from the WWC plasma through the cool material at the apex and/or the 
undisturbed primary wind, and the absorber which lies far outside the binary system. Iron K fluorescence 
emission is produced by scattering in the primary wind. One possibility for the CCE, where the CCE is 
produced in outer shock, is shown. Emission from the CCE only passes through the outer absorber beyond 
the binary system. 



companion's wind and thus reduce the emission 
from the colliding wind shock for a brief period. 
In optical observations, the He II A4687 emission 
line, which is believed to arise near the WWC 
shock, decreased in conce rt with the X-ray dro 
(|Steiner fc Daminelil Ex)! 



Martin et al. (200 



argued that the disappearance of this line in a di- 
rect view of the star, and in a re flected polar view 
(the so-called "FOS4" position, IStahl et al.ll2005l ) 
indicated a waning of the shocked gas during the 
minimum. 

Below, we examine some of the behaviors of the 
electron temperature, iron line profile, iron fluores- 
cence line variation and lack of variation at E > 
7 keV during the minimum from the 2003 X-ray 
spectra to try to decide between these two mech- 
anisms. 

6.3.1. Electron Temperature 

The electron temperature of the hottest mate- 
rial did not change significantly during the 2003 
minimum, as the normalized spectra at the bot- 
tom of Figure [7] clearly show. The constancy 
of the electron temperature of the hottest ma- 
terial argues against radiative braking, shock 



instability or accretion models which predict 
that the plasma temperature should decrease. 
lAkashi et alj ( 20061 ) 's model in particular predicts 
complete shutdown of the WWC activity during 
periastron passage, and therefore could only be 
consistent with a constant electron temperature 
if an alternative X-ray emission mechanism dur- 
ing the minimum could produce plasma at nearly 
the same temperature as the WWC emission. If 
the shock actually disappears during the mini- 
mum, then the constancy of the high-energy X- 
ray continuum requires that a small part of the 
3000 km s _1 secondary wind continues to pro- 
duce a shock by encountering material directly at 
normal incidence. 

6.3.2. The Iron Line Profile 

The variable excess emission on the red side 
of the Fe XXV line can be explained by non- 
equilibrium ionization (NEI) effects in the high- 
temperature plasma, and variations in the strength 
of this feature means that the ion temperature 
of some of the high-temperature plasma in the 
colliding wind shock changed during the mini- 
mum. This might favor mechanisms in which the 
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intrinsic emission from the shock declines with 
time near periastron, though the NEI effect itself 
does not reduce emissivity of the continuum. The 
observed excess at ^6.55 keV is consistent with 
emission lines of around Fe XXI, present at an 
ionization time scale r = nt ~ 10 10 ~ n s cm~ 3 
in the ~3.3 keV plasma, where n is the gas den- 
sity and t is the time after the gas starts heating. 
Since the gas density n oc D~ 2 , it should be larger 
by more than two orders of magnitude at perias- 
tr on than at apastron (u sing the orbital elements 
of I Corcoran et aL 2001a[ ). Thus near periastron 
the plasma should re ach ionizatio n equilibrium 
much faster (~15 min, Soker 2005) than at apas- 
tron. Perhaps near periastron, both the cooling 
time-scale and the equilibrium time-scale decrease 
relative to their values at apastron, but the cooling 
time-scale decreases faster. This means that near 
periastron the fraction of non-equilibrium gas to 
gas in equilibrium is larger than it is near apastron. 
However, since the radiative cooling time-scale has 
the same dependence on gas density as the ioniza- 
tion equilibrium time-scale (oc n -1 ), an additional 
cooling mechanism (such as adiabatic expansion 
or conduction) needs to be effective near perias- 
tron in order to decrease the cooling time-scale 
faster than the equilibrium time-scale. 

6.3.3. Iron Fluorescence Line Variations 

As noted in the previous section, the EWs of 
the Fe K fluorescence line stayed below ~200 eV 
during the 2003 minimum, and therefore the ob- 
served X-ray emission is direct emission, with only 
a small scattered component. The Fe fluorescence 
line flux decreased by a factor of ~50 during the 
2003 minimum (Figure [T3|) . This could be pro- 
duced by reducing the number of E >7.1 keV X- 
ray photons from the source seen by the fluorescing 
region, or by reducing the number of photons from 
the fluorescing region which reach the observer. If 
we assume that the fluorescing region extends far 
from the stars then during the 2003 minimum ei- 
ther the X-ray source is nearly entirely covered by 
obscuring material along the line of sight to the 
fluorescing region, or the intrinsic X-ray emissivity 
has declined so that fewer ionizing photons reach 
the fluorescing medium. This picture is consistent 
with a similar variation of He II A4687 emission be- 
tween the direct emission from the central WWC 
region and polar emission scattered from the SE 



lobe (|Stahl et al.ll2005HMartin et alj|200§) . 



6.3.4- The Lack of Variation at E > 7 keV Dur- 
ing the Minimum 

The hard band X-ray emission above ~7 keV 
reached its lowest value after the 2003 minimum 
started, recovered by a factor of 3 after XMM 030 722 
and did not vary for at least ^3 weeks from 
XMM 030802 to CXO 30828- Since RXTE light 
curves showed a similar trend in the previous cy- 
cle, the time variation seems to be phase repeat- 
able. The observed X-ray flux after XMM030722 
was ~5% of the flux at XMM 306i3- It is very 
difficult for an eclipse of the X-ray emitting re- 
gio n by an opaque disk (similar to that proposed 
by llshibashi et al. 19991 for example) , to keep the 
observed flux constant for so long a time, since 
the stellar motion would be so fast near perias- 
tron that the emitting region moves beyond the oc- 
culter on ~day timescales, causing a short eclipse 
and a rapid rise in X-ray brightness. Perhaps a 
"leaky absorber" eclipse, in which the absorbing 
medium is very extended and composed of opti- 
cally thick clumps, could better describe the be- 
havior of the observed flux during the minimum. 
Models in which the shock cools may not be con- 
sistent with the constancy of the E > 7 keV emis- 
sion either, unless some residual emission at high 
energies could be generated after the shock cools. 

6.4. The Origin of the CCE Component 

The CCE component did not vary within the 
two months during the 2003 minimum when it 
could be observed. It apparently had the same 
spectrum during the 1998 minimum (as shown by 
our analysis of ASCA971224), and did not show any 
enhanced emission during the minimum in 1992. 
This suggests that the X-ray source is stable on 
time-scales of 5—10 years. The Chandra image 
during the 2003 minimum restricts the size of the 
emitting plasma to < 1", constraining the maxi- 
mum size of the CCE at <2000AU. The column 
density to this component, TVh ^5x10 22 cm~ 2 , 
is smaller than the column density to the vari- 
able emission component around the 2003 mini- 
mum, but each component shows nearly the same 
column density in January 2003 (before the 2003 
minimum) and near apastron. After correcting for 
absorption, the X-ray luminosity of the CCE is 
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~10 3 ergs s . This suggests that a stable lu- 
minous source is present outside of the absorbing 
material which lies in front of the colliding-wind 
source at periastron. 

One possibility is that the CCE component 
is produced within the wind of the companion 
star and/or r\ Car due to the inherent insta- 
bilities of their radiatively-driven winds. The 
X-ray properties of the CCE component - kT 
~1 kcV with little apparent X-ray variation - 
resemble t he X-ray properties of normal OB 
stars (e.g. Berghofer et al. Il997t ICorcoran et al 
1993h. For OB stars log L x /i bol w -6 to 



(jBerghofer et al. 1997 ). so that the X-ray lumi- 
nosity of the CCE corresponds to a bolomet- 
ric luminosity Lboi ~ iq40-42 er g S g -i^ This 
is somewhat larger th an the luminosity of the 
companion according to Verner et al. ( 2005h (Lboi 



but nea r the luminosity o f 
HnTier et al.ll200lh . 



~ 4 x 10 39 ergs s x ) 
r] Car (L bo i ~ 2 x 10 40 ergs s 

Alternatively, the CCE emission could be pro- 
duced in ambient gas shocked by a strong out- 
flow from the stellar winds (see Figure [T4|) . The 
maximum size of the CCE (r < 2000AU) means 
that it is inside the Little Homunculus nebula 
( Ishibashi et al.l I2003T ) . so that the CCE compo- 



nent might be produced by the collision of the 
polar wind from rj Car with the Little Homuncu- 
lus. The velocity of the polar wind from r/ Car 
is believed to be v ^1000 km s _1 ( Smith et al 



2003al ). which is sufficient to heat gas to kT 
~1 keV. If all the kinetic energy of this out- 
flow is converted to thermal energy, the required 
mass loss rate would have to be only ^3xl0 -8 
Mq yr _1 in order to match the observed emis- 
sion measure (~7.0xl0 56 cm -3 ), which is only 
3 x 10~ 4 — 3 x 10~ 5 of the mass loss rate of 
n Car (~1CT 3 -10- 4 tf w yr" 1 , Iffillier et alJl200ll 
Pittard fc Corcoran] 2002 ). The gas density inside 
the shell structure of the Little Homunculus should 
be less than the den sity of the shell, wh ich is 
10 6 < n H <10 7 cm" 3 (jlshibashi et al.ll2003l ). As- 
suming that radiative cooling dominates, the cool- 
ing timescale of the plasma is r ^3nfcT/n 2 A(T), 
where A(T) is the emissivity. With a plasma den- 
sity of n <10 7 cm -3 , and plasma temperature kT 
-1 keV, and A(lkeV) - 3.1 xl0~ 23 ergs s" 1 cm 3 , 
then r >0.5 years, consistent with our result. 

A third alternative is that the CCE might be 
associated with the outflowing material from the 



WWC that collides with cold circumstellar gas 
and is reheated. The required mass loss rate 
streaming out through the WWC is given by M = 
2L x /eV 2 ~10~ 8 /(e x (F/2000km s" 1 ) 2 ) M Q yr" 1 , 
where e is the efficiency to convert the kinetic en- 
ergy to X-ray emission, and V is the velocity of 
the outflow, which should be >1,000 km s _1 from 
the CCE temperature. Since the most of the com- 
panion's wind goes into the shock cone, the M of 
the outflowing material sho uld be ~10 -5 M Q yr^ 1 
(jPittard fc Corcoran! 120021 ). and e should only be 
-0.1%. 

7. Conclusions &: Summary 

We have observed r\ Car at key phases around 
the X-ray minimum in 2003 using the XMM- 
Newton and Chandra satellites and studied the 
variation of the X-ray spectral shape. We also 
compared the behavior of the X-ray emission dur- 
ing the 2003 minimum with observations during 
earlier minima made by ASCA, BeppoSAX and 
ROSAT. The results are summarized as follows. 

1. During the minimum, the hard and medium 
band emission from the central point source 
decreased to 1—5% of the flux outside the 
minimum, but the emission did not disap- 
pear. The X-ray flux varied in an appar- 
ently stepwise way near the middle of the 
minimum, suggesting that the X-ray mini- 
mum has two states. The light curve within 
any single observation did not show signifi- 
cant variability on timescales of a few ksec, 
supporting that the hard X-ray emission is 
produced in a WWC. 

2. The spectrum of the variable emission com- 
ponent changed markedly throughout the 
X-ray minimum. The medium band emis- 
sion decreased gradually relative to the 
hard emission. Assuming an absorbed IT 
model for the entire spectrum, the absorp- 
tion to the central source increased from 
5xl0 22 cm~ 2 to 4xl0 23 cm -2 . The hard 
band flux went up prior to the X-ray min- 
imum, and decreased abruptly during the 
2003 minimum, without changing the spec- 
tral slope significantly, suggesting that the 
plasma electron temperature is constant 
around 3 keV. 
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The line flux ratios between the helium- 
like ions and hydrogen-like ions of Si, S, 
and Fe seen in the spectra before the 2003 
minimum indicate a similar type of multi- 
te mperature plasma as se en near apastron 
bv lCorcoran et alj (|2001bl ). 



4. The values of JVh derived from the spectra 
above 5 keV (JVh[>5 keV]) are larger than 
those values of Nn derived from the entire 
spectra (TVh [whole]) or from the spectra be- 
low 5 keV (JVh[<5 keV]), but show less vari- 
ation than N-a [whole] or iV H [<5 keV]. This 
suggests that the hotter plasma has stronger 
absorption, i.e. emission from r\ Car is multi- 
temperature, and each emission component 
is associated with an independent absorption 
component. This means that the hotter and 
cooler emission originate in spatially distinct 
regions. The hardest emission, which is be- 
lieved to arise from material near the stag- 
nation point of the WWC region, perhaps 
shows the effects of a build-up of material as 
the colliding wind shock plows through the 
primary wind. 

5. The Fe XXV line profile showed a small 
blue excess, and a strong red excess before 
the 2003 minimum. The blue excess may 
be due to Doppler shifting of the line cen- 
troid or line broadening, and if so, the ob- 
served velocity is ~4000 km s _1 . The red 
excess may indicate that some of the plasma 
is no longer in ionization equilibrium. The 
line showed a stronger red excess during the 
2003 minimum than outside the minimum, 
which we interpret as decrease in the cool- 
ing time scale of the plasma near periastron. 
Interestingly, we note that an outburst spec- 
trum of the dwarf nova SS Cyg taken with 
the Suzaku satellite showed a similar red ex- 
cess on the Fe K line profile accompanied 
with a decrease in the hard band X-ray flux 
(jlshida et al.ll2007h . 



6. We find evidence for a hard tail at E > 9 
keV, at about an order of magnitude larger 
than the background level, in spectra out- 
side the X-ray minimum. This either in- 
dicates the presence of very hot plasma in 
the system or non-thermal emission perhaps 
produced by first-order Fermi acceleration of 



particles which then Compton-upscatter UV 
seed photons from the stellar photospheres 
and/or X-ray photons from the shock. 

7. The values of Nh and kT of the X-ray 
spectrum at corresponding phases in pre- 
vious cycles are very similar to those during 
the 2003 cycle, except during the recovery 
phase. In particular, the ASCA spectrum 
during the 1997/98 minimum (ASCA971224) 
is almost identical to the spectrum from 
CXOXMM 030 7. The difference in JV H during 
the recovery, as measured by Chandra and 
BeppoSAX could be produced by an increase 
in the density of the wind from 77 Car. 

8. The EWs of the Fe fluorescence line at 
6.4 keV are always below the values ex- 
pected from ./Vh[>5 keV] assuming a spher- 
ically symmetric absorber. This is probably 
because the absorber is not spherically sym- 
metric and is consistent with the lack of 
scattering matter in the lower-density wind 
of the companion if the full opening angle of 
the shock cone i s ~ 120°, consistent with th e 
value derived bv lPittard fe Corcoran ( 2002 ). 



9. The X-ray minimum is caused by the appar- 
ent decrease of the plasma E.M. An eclipse 
model seems to explain the observed con- 
stancy of the plasma electron temperature, 
while the small EWs of the Fe fluorescence 
line and the deformation of the Fe XXV He- 
like triplet might indicate a real decline in 
the amount of hot gas in the colliding wind 
shock. It might also indicate that nearly all 
of the shock is hidden by an optically thick 
but slightly porous stellar wind. 

10. We discovered a constant emission (CCE) 
component very near the central X-ray 
source. The CCE component can be seen 
between 1—3 keV during the X-ray mini- 
mum; at other times the emission from the 
colliding wind source overwhelms it. This 
component showed no variation on either 
short timescales within any observation or 
between observations for about two months 
during the 2003 minimum, and was appar- 
ently the same during the 1997/98 and 1992 
minima. The CCE spectrum exhibited emis- 
sion lines of Si XIII and S XIV, and can 
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be fit by an absorbed IT model with Nn 
~5xl0 22 cm~ 2 , kT ~1.1 keV and absorp- 
tion corrected Lx ~10 34 ergs s _1 . Chandra 
data suggest that the emission is point-like, 
restricting the size of the source <2000AU 
at 2.3 kpc. This emission might be produced 
by collisionally-heated shocks from the fast 
winds from rj Car or perhaps the fast moving 
outflow from the wind-wind collision with 
the ambient gas, or might be associated with 
shocked gas which is intrinsic to the wind of 
T) Car. 

11. X-ray emission within 3' from rj Car can be 
reproduced by 2 emission components with 
interstellar absorption: cool plasma with kT 
^0.25 keV (perhaps coming from the diffuse 
emission and HDE 303308), and hot plasma 
with fcT~2.6 keV or power-law T ~ 2.5, per- 
haps coming from a cluster of young stars. 
Correcting for these emission components is 
important in understanding the X-ray spec- 
trum of rj Car derived from lower spatial res- 
olution X-ray observatories like ASCA and 
Suzaku. 

12. Emission from the outer ejecta did not vary 
significantly on timescales of several years. 
The slope of the highest energy emission 
associated with the spectrum of the outer 
ejecta can be reproduced with kT ~0.6 keV, 
which is consistent with the highest temper- 
ature derived from analysis of the em ission 
line spectrum (jLeutenegger et aL 20031 ). 



Currently, there are few reasons to suspect that 
the variable X-ray emission is not driven by shocks 
produced by the WWC from binary stars. For the 
first time, the spatially-resolved observations of 
i] Car obtained by XMM-Newton and Chandra has 
revealed the details of the plasma and absorption 
components which are directly associated with the 
W WC plasma. While the ear lier spectral analysis 
of IPittard fe Corcoran! (|2002h did not account for 
the different emission components identified in this 
work, their fit was largely determined by emission 
above 2 keV which we have shown is dominated 
by emission from the WWC at the phase of their 
analysis. 

Still the biggest mystery is the cause of the X- 
ray minimum, which is difficult to fully explain by 



occultation by a smooth wind from the primary 
star. If the X-ray activity really switches off dur- 
ing the minimum, a fading of the X-ray emission 
can be imaged in principle using the scattered X- 
ray emission from the Homunculus Nebula. The 
intensity of the Fe K-shell fluorescence line just af- 
ter the onset of the 2003 minimum should change 
differently between the eclipsing model and the 
fading model, and perhaps detailed monitoring of 
this line by spatially resolved X-ray spectral obser- 
vations could decide the issue. On the other hand, 
deformation of the Fe XXV line profiles is unusual 
in the context of stellar systems, though similar 
deformation has been noted in iron lines arising 
from regions of strong gravity near black holes. 
This feature requires further study with high res- 
olution spectrometers such as an X-ray micro- 
calorimeter to produce sufficient spectral resolu- 
tion at 7 keV. 
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Misaki and K. Mukai for useful comments. We 
greatly appreciate the XMM-Newton and Chan- 
dra review committee's positive response for the 
TOO observation proposal. This work is per- 
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NASA/GSFC, and supported by Chandra US 
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